Measurements of Direct CP Violation, CPT Symmetry, 
and Other Parameters in the Neutral Kaon System 
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We present a series of measurements based on Kl,s —> it + tt~ and Kl,s — * 7r°7i"° decays collected 
in 1996-1997 by the KTeV experiment (E832) at Fermilab. We compare these four K — > nn decay 
rates to measure the direct CP violation parameter Re(e'/e) = (20.7 ± 2.8) x 10~ 4 . We also test 
CPT symmetry by measuring the relative phase between the CP violating and CP conserving 

decay amplitudes for K — > tt + it~ (<^>_| ) and for K — » 7r°7r° ((poo)- We find the difference between 

the relative phases to be A<fi = <j>oo — 4>+- = (+0.39 ± 0.50)°, and the deviation of from 

the superweak phase to be (j>-\ — (f>sw = (+0.61 ± 1.19)°; both results are consistent with CPT 

symmetry. In addition, we present new measurements of the Kl-Ks mass difference and Ks 
lifetime: Am = (5261 ± 15) x 10 6 hs' 1 and r S = (89.65 ± 0.07) x 10~ 12 s . 

PACS numbers: 11.30.Er, 13.25.Es, 14.40.Aq 



I. INTRODUCTION 

The discovery of the Kt, — > tt + tt~ decay in 1964 Q| 
demonstrated that CP symmetry is violated in weak in- 
teractions. Subsequent experiments showed that the ef- 
fect is mostly due to a small asymmetry between the 
K° — » K° and K° — > K° transition rates, which is re- 
ferred to as indirect CP violation. Over the last three 
decades, significant effort has been devoted to searching 
for direct CP violation in a decay amplitude. 

Direct CP violation can be detected by comparing the 
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level of CP violation for different decay modes. The pa- 
rameters e and d are related to the ratio of CP violating 
to CP conserving decay amplitudes for K — » 7r + 7r~ and 
K -» 7r°7r°: 
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e is a measure of indirect CP violation, which is common 
to all decay modes. If CPT symmetry holds, the phase 
of e is equal to the "superweak" phase: 

<j) SW = tan" 1 (2Am/Ar) , (2) 

where Am = ttil — ms is the K^-Ks mass difference and 
Ar = Ts — r_L is the difference in the decay widths. 
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The quantity e' is a measure of direct CP violation, 
which contributes differently to the 7r + 7r~ and 7r°7r° de- 
cay modes, and is proportional to the difference be- 
tween the decay amplitudes for K° — > tt + it~ (ir°ir°) and 
if — ► 7T + 7r _ (7r 7r°). Measurements of irir phase shifts || 
show that, in the absence of CPT violation, the phase of 
e' is approximately equal to that of e. Therefore, i?e(e'/e) 
is a measure of direct CP violation and Im(e'/ e) is a mea- 
sure of CPT violation. 

Experimentally, Re(e'/e) is determined from the double 
ratio of the two pion decay rates of Kl and ifg: 

T(K L -> 7T+7T-) /T(K S -> 7T+7T-) 
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For small |e'/e|, Im(e'/e) is related to the phases of 77.) 
and 7700 by 
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The Standard Model accommodates both direct and 
indirect CP violation ||, [|. Unfortunately, there are 
large hadronic uncertainties in the calculation of Re(e'/e). 
Most recent Standard Modelpredictions for Re(e'Je) are 
less than 30 x 10~ 4 § R | g 00 H El II- 
The superweak model [n7] , proposed shortly after the dis- 
covery of Kl — > 7r + 7r _ , also accommodates indirect CP 
violation, but not direct CP violation. 

Previous measurements have established that Re(e'/e) 
is non-zero |l8l [l9|, |2l[]. This paper reports an im- 
proved measurement of Re(e'/e) by the KTeV Experi- 
ment (E832) at Fermilab. This measurement is based 
on 40 million reconstructed if — > irir decays collected in 
1996 and 1997, and represents half of the total KTeV data 
sample. The 1996+1997 sample is four times lar ger than, 
and contains, our previously published sample |20j. We 
also present measurements of the kaon parameters Am 
and ts, and tests of CPT symmetry based on measure- 
ments of A(f> and <f>^ — 4>sw- 

The outline of the paper is as follows. Section [H] de- 
scribes the KTeV measurement technique, incl udin g de- 
tails about the neutral hadron beams (Sec. [II B| ) and 
the detector used to identify if — ► nn decays (Sec. II C). 



The detector description also includes the calibration 
procedures and the performance characteristics relevant 
to the Re(e'/e) measurement. The Monte Carlo simu- 
lation of the kaon beams and detector is described in 



Section II E . Section III explains the reconstruction and 
event selection for the if — > 7r + 7r~ and if — > 7r°7r° de- 
cay modes, and also the background subtraction for each 
mode. Section IV describes the acceptance correction, 
and the fit used to extract i?e(e'/e) and other physics 
parameters. Each of these sections is followed by a dis- 
cussion of systematic uncertainties related to that part 
of the Re(e'/e) measurement. Section [v] presents the 
Re(e'/e) result along with several crosschecks. Finally, 



Section VI presents our measurements of the kaon pa- 
Additional details of 



rameters Am, ts, A</>, and 
the work presented here are given in 



II. MEASUREMENT TECHNIQUE AND 
APPARATUS 



A. Overview 

The measurement of Re(e'/e) requires a source of Kl 
and Ks decays, and a detector to reconstruct the charged 
(7r + 7r~) and neutral (7r 7r°) final states. The strategy 
of the KTeV experiment is to produce two identical Kl 
beams, and then to pass one of the beams through a 
thick "regenerator." The beam that passes through the 
regenerator is called the regenerator beam, and the other 
beam is called the vacuum beam. The regenerator creates 
a coherent \Kl) + p\Ks) state, where p is the regener- 
ation amplitude chosen such that most of the K — > tttt 
decay rate downstream of the regenerator is from the 
Ks component. The measured quantities are the num- 
bers of if — > 7r + 7r~ and K — > 7r°7r° decays in the vac- 
uum and regenerator beams. The vacuum-to-regenerator 
"single ratios" for if — * 7r + 7r~ and K — > 7r°7r° decays are 

proportional to \r]-\ / p\ 2 and |?7oo/p| 2 , and the ratio of 

these two quantities gives Re(e'/e) via Eq. || (also see 
Appendix |Aj) . The effect of Ks-Kl interference in the 
regenerator beam is accounted for in a fitting program 
used to extract Re(e'/e). 

To reduce systematic uncertainties related to left-right 
asymmetries in the detector and beamline, the regenera- 
tor position alternates between the two beams once per 
minute. Decays in both beams are collected simultane- 
ously to reduce sensitivity to time-dependent variations 
in the beam intensity and in detector efficiencies. The 
fixed geometry of the beamline elements, combined with 
alternating the regenerator position, ensures a constant 
vacuum-to-regenerator kaon flux ratio. 

K — ► 7r + 7r~ decays are detected in a spectrometer con- 
sisting of four drift chambers and a dipole magnet; the 
well-known kaon mass is used to determine the momen- 
tum scale. The four photons from K — > 7r°7r° decays are 
detected in a pure Cesium Iodide (Csl) calorimeter; elec- 
trons from Kl — > Tr ± e T i' decays (if e 3) are used to cal- 
ibrate the Csl energy scale. An extensive veto system 
is used to reject events coming from interactions in the 
regenerator, and to reduce backgrounds from kaon de- 
cays into non-7T7r final states such as Kl — ► ^p^v and 
K L tt'W. 

A Monte Carlo simulation is used to correct for the ac- 
ceptance difference between K — > tttt decays in the two 
beams, which results from the very different Kl and if 5 
lifetimes. The simulation includes details of detector ge- 
ometry and efficiency, as well as the effects of "acciden- 
tal" activity from the high flux of particles hitting the 
detector. The decay-vertex distributions provide a crit- 
ical check of the simulation. To study the detector per- 
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FIG. 1: Top view of the KTeV beamline, showing the BeO 
target, absorbers, collimating system, Mask Anti, and regen- 
erator. The common absorber is a 7.62 cm thick Pb slab 
followed by a 52.3 cm slab of beryllium. The movable ab- 
sorber is a 45.7 cm slab of beryllium that alternates beam 
positions once every minute along with the regenerator. The 
crossover absorber is a 2 m thick iron slab. Note the different 
horizontal and vertical scales. 



formance, and to verify the accuracy of the Monte Carlo 
simulation, we also collect samples of Kl — > 7r ± e =F ^ and 
Kl — » tt°tt°tt decays with much higher statistics than 
the K — > 7T7T signal samples. 



B. The Kaon Beams 

Neutral kaons are produced by a proton beam hitting 
a fixed target (Fig. |l). The Fermilab Tevatron provides 
3 x 10 12 800 GeV/c protons in a 20 s extraction cycle 
("spill") once every minute. The proton beam has a 
53 MHz RF micro-structure such that protons arrive in 
~ 1 ns wide "buckets" and at 19 ns intervals. The bucket- 
to-bucket variations in beam intensity are typically 10%. 

The target is a narrow beryllium oxide (BeO) rod, 3 x 
3 mm 2 in the dimensions transverse to the beam, and 
30 cm long corresponding to about one proton interaction 
length. The proton beam is incident on the target at 
a downward angle of 4.8 mrad with respect to the line 
joining the target and the center of the detector. This 
targeting angle is chosen as a compromise between higher 
kaon flux at small angles, and a smaller neutron-to-kaon 
ratio (n/K) at large angles. 

The center of the BeO target defines the origin of the 
KTeV right-handed coordinate system. The positive z- 
axis is directed from the target to the center of the de- 



tector. The positive y-axis is directed up. 

The particles produced in the BeO target include very 
few kaons compared to other hadrons and photons. A 
~ 100 m long beamline is used to remove unwanted par- 
ticles and to collimate two well-defined kaon beams. In 
this beamline, charged particles are removed with sweep- 
ing magnets, photons are absorbed by a 7.62 cm Pb slab 
located at z = 19 m, and most of the hyperons decay near 
the target. To reduce the neutron-to-kaon ratio, neutrons 
(kaons) are attenuated by a factor of 4.6 (2.6) in a beryl- 
lium absorber that is common to both beams. An extra 
"movable" absorber, synchronized with the regenerator 
position, provides an additional attenuation factor of 3.8 
(2.3) for neutrons (kaons) in the regenerator beam. 

Each neutral kaon beam is defined by two collimators: 
a 1.5 m long primary collimator [z = 20 m), and a 3 m 
long "defining" collimator (z = 85-88 m) that defines 
the size and solid angle of each beam. Each collimator 
has two square holes, which are tapered to reduce scat- 
tering. A "crossover" absorber at z = 40 m prevents 
kaons that scatter in the upstream absorbers from cross- 
ing over into the other beam. At the defining collimator, 
the two beams have the same size (4.4 x 4.4 cm 2 ) and 
solid angle (0.24 /istr). The beam centers are separated 
by 14.2 cm, and the horizontal angle between the two 
beams is 1.6 mrad. 

The two beams pass through an evacuated volume, 
held at 10 -6 Torr, which extends from 28 m to 159 m 
from the target. This evacuated region includes the 
110 m to 158 m range used in the analysis. The down- 
stream end of the evacuated volume is sealed with a 
0.14% radiation length (Xq) vacuum-window made of 
kevlar and mylar. 

Most of the Ks component decays near the BeO target. 
Downstream of the defining collimator, the small "target- 
Ks" component that remains in the vacuum beam in- 
creases the K — > 7T7T decay rate by 0.4% compared with 
a pure Kl beam of equal intensity. This contribution 
from target-i^g is essentially zero for kaon momenta be- 
low 100 GeV/c and rises to 15% at 160 GeV/c. 

In the regenerator beam, most of the K — > irir decays 
are from the regenerated Ks component. Decays from 
the Kl component and from K$-Kl interference account 
for about 20% of the K — > 7T7t decay rate, ranging from 
5% near the regenerator to 90% at the vacuum-window. 

The composition of the neutral hadron beams is as fol- 
lows. At 90 m from the BeO target, the vacuum beam has 
a 2.0 MHz flux of kaons with n/K — 1.3, and an average 
kaon momentum of about 70 GeV/c. Upstream of the 
regenerator, the kaon flux is 0.9 MHz with n/K = 0.8; 
downstream of the regenerator the flux of unscattered 
kaons is 0.15 Mhz, or 13 times smaller than in the vac- 
uum beam. The flux of hyperons at 90 m from the tar- 
get is about 1 kHz in the vacuum beam. In addition to 
hadrons, there are muons that come from the BeO tar- 
get, the proton beam dump, and kaon decays; the total 
muon flux hitting the KTeV detector is about 200 kHz. 
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FIG. 2: Top view of the KTeV (E832) detector. The evacuated decay volume extends to z = 159 m. 



C. The Detector 

Kaon decays downstream of the defining collimator 
(z = 85 m, Fig. are reconstructed by the KTeV de- 
tector (Fig. ||), which includes a magnetic spectrometer, 
Csl calorimeter, and veto system. Downstream of the 
vacuum- window, the space between detector components 
is filled with helium to reduce interactions from the neu- 
tral beam, and to reduce multiple scattering and photon 
conversions of decay products. The total amount of ma- 
terial upstream of the Csl calorimeter corresponds to 4% 
of a radiation length; about 60% of the material is in the 
trigger hodoscope just upstream of the Csl calorimeter, 
and 10% of the material is upstream of the first drift 
chamber. Each of the two neutral beams passes through 
holes in the Mask Anti veto, trigger hodoscope, and Csl 
calorimeter. The beams finally strike a beam-hole pho- 
ton veto 5 meters downstream of the Csl. The following 
sections describe the detector components in more detail. 



1. Spectrometer 

The KTeV spectrometer includes four drift chambers 
(DCs) that measure charged particle positions. The two 
downstream chambers are separated from the upstream 
chambers by a 3 x 2 m 2 aperture dipole magnet. The 
magnet produces a field that is uniform to better than 1% 
and imparts a 0.41 GeV/c kick in the horizontal plane. 
During data taking, the magnet polarity was reversed 
every 1-2 days. 

The DC planes have a hexagonal cell geometry formed 



by six field-shaping wires surrounding each sense wire 
(Fig. [|). The cells are 6.35 mm wide, and the drift ve- 
locity is about 50/mi/ns in the 50-50 argon-ethane gas 
mixture. The maximum drift time across a cell is 150 ns, 
and defines the width of the "in-time" window. A cham- 
ber consists of two planes of horizontal wires to measure y 
hit positions, and two planes of vertical wires to measure 
x hit positions; the two ir-planes, as well as the two y- 
planes, are offset by one half-cell to resolve the left-right 
ambiguity. The x and y hits cannot be associated to each 
other using only drift chamber information, but c an be 
matched using Csl clusters as explained in Section III B . 

There are a total of 16 planes and 1972 sense wires in 
the four DCs. The transverse chamber size increases with 
distance from the target. The smallest chamber (DC1) 
is 1.26 x 1.26 m 2 and has 101 sense wires per plane; the 
largest chamber (DC4) is 1.77 x 1.77 m 2 and has 140 
sense wires per plane. 

Lecroy 3373 multi-hit Timc-to-Digital Converters 
(TDCs) are used to measure the drift times relative to the 
RF-synchronized Level 1 trigger. The TDC resolution is 
0.25 ns, which contributes 13 fim to the position reso- 
lution. The total TDC time window is 2.5 times longer 
than the in-time window, and is centered on the in-time 
window. The track reconstruction software uses only the 
earliest in-time hit on each wire. Hits prior to the in-time 
window are recorded to study their influence on in-time 
hits. 

Each measured drift time, t, is converted into a drift 
distance, x, with a non-linear x(t) map. The maps are 
measured separately for each of the sixteen planes using 
the uniform hit-illumination across each cell. The x(t) 
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FIG. 4: Deviation of the sum of drift distances (SOD) 
from the nominal cell size of 6.35 mm, for the y-views us- 
ing a sample of K — » -k + -k~ decays. The arrows indicate the 
±1 mm requirement for a good hit-pair. The data are shown 
as dots. The histogram is the Monte Carlo (MC) prediction 
(Sec. Ee|). 



FIG. 3: Drift chamber hexagonal cell geometry showing six 
field wires (open circles) around each sense wire (solid dots). 
The vertical dashed lines denote the boundaries of the "offset- 



cells" used in the track-separation cut (see Section 1IIB). For 
example, two offset-cells are between tracks 1 and 2. 



calibrations are performed in roughly 1-day time periods. 

A charged-particle track produces a hit in each of the 
two x and y planes. The two x hits (or two y hits) in 
each plane are referred to as a "hit-pair." For a track that 
is perpendicular to a drift chamber with perfect resolu- 
tion, the sum of drift distances (SOD) from each hit-pair 
would equal the cell width of 6.35 mm. The measured 
SOD distribution is shown in Fig. |J For inclined tracks, 
an angular correction is applied to the SOD. The track- 
finding software requires the SOD to be within 1 mm 
of the cell width. The mean SOD is stable to within 
10 fim during the run. The single-hit position resolution 
is typically 110 ^m, corresponding to a SOD resolution 
of 150 /im and a hit-pair resolution of 80 fim. Using the 
tracking algorithm described in Section IIIB, the mo- 
mentum resolution is <r p /p ~ [1.7© (p/14)] x 10~ 3 , where 
p is the track momentum in GeV/c. 

The average inefficiency for reconstructing a hit-pair 
is 3.7%. Delta rays and accidental hits contribute to 
the low-side tail of the SOD distribution (Fig. [|), and 
result in a hit-pair loss of 0.5% and 0.7%, respectively. 
The remaining 2.5% loss is from missing hits and from 
SOD values more than 1 mm larger than the cell size 
(high-side tail in Fig. ^) . More details on sources of hit- 
pair inefficiency are given in the description of the Monte 
Carlo simulation (Sec. [IE 2). 



be known. This alignment is determined in roughly 1-day 
periods using data samples described below. 

The transverse alignment of the drift chambers is based 
on muons from dedicated runs with the analysis magnet 
turned off. The muon intensity is raised to ~ 1 MHz by 
reducing the field in the upstream sweeping magnets, and 
the neutral hadron beam is absorbed by a 2 m long steel 
block placed in the beam at z = 90 m. These muon runs 
were performed every 1-2 days when the magnet polarity 
was reversed. Software calibration results in a transverse 
alignment of each x and y plane to within 10 /im, and 
relative rotations known to within 20 /irad. 

The transverse target position is determined with a 
precision of 35 ^m using K — » tt + tt~ decays in the vac- 
uum beam and projecting the reconstructed kaon trajec- 
tory back to the target. The Csl offset relative to the drift 
chambers is measured using electrons from Kl 
decays. 



To reconstruct the trajectories of charged particles ac- 
curately, the alignment of the drift chambers relative to 
each other, to the target, and to the Csl calorimeter must 



2. The Csl Calorimeter 

The KTcV electromagnetic calorimeter consists of 3100 
pure Csl crystals as shown in Fig. |[ There are 2232 2.5 x 
2.5 cm 2 crystals in the central region, each viewed by a 
1.9 cm Hamamatsu R5364 photomultiplier tube (PMT). 
There are also 868 5 x 5 cm 2 crystals, each viewed by a 
3.8 cm Hamamatsu R5330 PMT. The transverse size of 
the calorimeter is 1.9 x 1.9 m 2 , and the length of each 
crystal is 50 cm (27 Xq). Two 15 x 15 cm 2 carbon fiber 
beam pipes allow the few MHz of beam particles to pass 
through the calorimeter without striking any material. 

The crystals were individually wrapped and tested 
with a 137 Cs source to ensure that the response over the 
length of each crystal is uniform to within ~ 5% [ p3| . 
This longitudinal uniformity requirement is necessary to 
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< 1.9 m > 

FIG. 5: Beamline view of the KTeV Csl calorimeter, show- 
ing the 868 larger outer crystals and the 2232 smaller inner 
crystals. Each beam hole size is 15 x 15 cm 2 and the two 
beam hole centers are separated by 0.3 m. 



obtain sub-percent energy resolution for electrons and 
photons. 

The Csl scintillation light has two components: (i) a 
fast component with decay times of 10 ns and 36 ns, 
and maximum light output at a wavelength of 315 nm; 
(ii) a slow component with a decay time of ~ 1/is and 
maximum light output at 480 nm. To reduce accidental 
pile-up effects from the slow component, a Schott UV fil- 
ter is placed between the crystal and the PMT. The filter 
reduces the total light output by ~ 20%, but increases 
the fast component fraction from about 80% to 90%. The 
average light yield with the filter gives 20 photo-electrons 
per MeV of energy deposit. The Csl signal components 
discussed above include the effect of the PMT spectral 
response. 

Digitizing electronics are placed directly behind the 
PMTs to minimize electronic noise (< 1 MeV). Each 
PMT is equipped with a custom 8-range digitizer to in- 
tegrate the charge delivered by the PMT. This device, 
which is referred to as a "digital PMT" (DPMT) @, 
has 16 bits of dynamic range with 8-bit resolution, and 
allows the measurement of energies from a few MeV to 
100 GeV. In 1997, the digitization and readout oper- 
ated at the Tevatron RF frequency of 53 MHz, and the 
PMT signal integration time was 114 ns (6 RF "buck- 
ets"), which permitted collection of approximately 96% 
of the fast scintillation component. In 1996, the readout 
frequency was 18 MHz (RF/3) and the integration time 
was a factor of two longer than in 1997. 

To convert measurements of integrated charge to en- 



ergy, a laser system is first used to calibrate the response 
from each DPMT. Then momentum-analyzed electrons 
from Kl — > 7r ± e T ^ decays are used to calibrate the en- 
ergy scale of each channel. 

The laser system consists of a pulsed Nd:YAG laser, 
a fluorescent dye, and a variable filter to span the full 
dynamic range of the readout system. The laser light 
(360 nm) is distributed via 3100 quartz fibers to each in- 
dividual crystal; the light level for each quadrant of the 
calorimeter is monitored with a PIN diode read out by a 
20-bit ADC. Throughout the data taking, special hour- 
long laser scans to calibrate the readout system were per- 
formed roughly once per week during periods when beam 
was not available. Using these laser scan calibrations, 
deviations from a linear fit of the combined DPMT plus 
PMT response versus light level are less than 0.1% (rms) 
for each channel. During nominal data-taking, the laser 
operated at a fixed intensity with 1 Hz repetition rate 
to correct for short-term gain drifts that were typically 
< 0.2% per day. 

To determine the energy scale in the calorimeter, we 
collected 600 million Kl — ► 7r ± e =F ^ decays during the ex- 
periment; this number of events allows us to determine 
the energy scale for each channel with ~ 0.03% precision 
every 1-2 days. The electron energy is determined by 
summing energies from a 3 x 3 "cluster" of large crystals 
or a 7 x 7 cluster of small crystals centered on the crystal 
with maximum energy. The cluster energy is corrected 
for shower leakage outside the cluster region, leakage at 
the beam-holes and calorimeter edges, and for channels 
with energies below the ~ 4 MeV readout threshold. 

Figure shows the E/p distribution for electrons from 
K e 3 decays, where E/p is the ratio of cluster energy mea- 
sured in the Csl calorimeter to momentum measured in 
the spectrometer. To avoid pion shower leakage into the 
electron shower, the tt ± is required to be at least 50 cm 
away from the e T at the Csl. The E/p resolution has 
comparable contributions from both E and p. The Csl 
energy resolution is obtained by subtracting the momen- 
tum resolution from the E/p resolution, and is shown as 
a function of momentum in Fig. ^|b. The energy resolu- 
tion can be parameterized as oe/E ~ 2%/yE © 0.4%, 
where E is in GeV; the resolution is 1.3% at 3 GeV, the 
minimum cluster energy used in the K — > 7r°7r° analysis, 
and is better than 0.6% for energies above 20 GeV. The 
momentum dependence of the mean E/p (Fig. ^p) shows 
that the average energy nonlinearity is 0.5% between 3 
and 75 GeV. This energy nonlinearity is measured for 
each Csl channel and used as a correction. 

Electromagnetic cluster positions are determined from 
the fraction of energy in neighboring columns and rows. 
The conversion from energy fraction to position is done 
using a map based on the uniform photon illumination 
across each crystal for K — > 7r°7r data. The average po- 
sition resolution for electrons is 1.2 mm for clusters in 
small crystals, and 2.4 mm for large crystals. 
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FIG. 6: (a) E/p distribution based on electrons from 
Kl — * Jre^v decays, where E is the energy measured in the 
Csl calorimeter and p is the momentum measured in the spec- 
trometer; (b) electron energy resolution vs. momentum (dots) 
obtained by subtracting the tracking resolution, shown by the 
smooth curve, from the _E/p-resolution; (c) E/p-mea,n vs. mo- 
mentum before the linearity correction is applied. 



3. The Regenerator 
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FIG. 7: (a) Diagram of the regenerator and (b) magnified 
view of the downstream end of the regenerator. PMTs are not 
shown. The dimension of each element is given in the figure. 
The kaon beam enters from the left. Vertical arrows above 
the lead and scintillator show the effective edge locations for 
K — > 7r°7r° and K — > -k + tt~ decays; the width of each arrow 
indicates the uncertainty. 



KTeV uses an active regenerator as a source of K$ de- 
cays. It consists mainly of 84 10 x 10 x 2 cm 3 scintillator 
modules (Fig. |^a). Each module is viewed by two photo- 
multiplier tubes (PMTs), one from above and one from 
below. The downstream end of the regenerator (Fig. [7b) 
is a lead-scintillator sandwich, which is also viewed by 
two PMTs. This last module of the regenerator is used 
to define a sharp upstream edge for the kaon decay region 
in the regenerator beam. 

At the average kaon momentum of 70 GeV/c, the mag- 
nitude of the regeneration amplitude is \p\ ~ 0.03. The 
isoscalar carbon in the regenerator accounts for about 
95% of the regeneration amplitude, which simplifies the 
model used to describe p when extracting physics param- 
eters (Sec. [fVC]). 

One can distinguish three main processes that con- 
tribute to Ks regeneration. These are (i) "coherent" 
regeneration, which occurs in the forward direction, (ii) 
"diffractive" regeneration, in which target nuclei do not 
disintegrate but kaons scatter at finite angle, and (iii) 
"inelastic" regeneration, characterized by nuclear break- 
up and often by production of secondary particles. Only 
the decays of coherently regenerated kaons are used in 
the Re(e'/e) analysis. The other processes are treated as 
background. 



The 170 cm length of the regenerator corresponds to 
about two hadronic interaction lengths. This length max- 
imizes coherent regeneration and suppresses diffractive 
regeneration p5[ . The diffractive-to-coherent ratio is 0.09 
for K — > 7T7T decays downstream of the regenerator, and 
is reduced by kinematic cuts in the analysis. The corre- 
sponding inelastic-to-coherent ratio is about 100. Since 
inelastic regenerator interactions typically leave energy 
deposits of a few MeV to 100 McV from the recoil nu- 
clear fragments, this source of background is reduced us- 
ing the regenerator PMT signals; events with more than 
8 MeV in any scintillator element are rejected. Inelastic 
interactions with production of secondary particles are 
suppressed further by other elements of the veto system. 
After all veto requirements, the level of inelastic scatter- 
ing is reduced by a factor of a few thousand, making its 
contribution smaller than that of diffractive scattering. 

The downstream edge of the regenerator defines the 
beginning of the regenerator-beam decay region for both 
the K — > 7r + 7r~ and K — ► 7r°7r° modes. A small fraction 
of decays inside the regenerator enters the signal sample 
because photons can pass through scintillator and lead, 
and because charged pions can exit the last regenerator 
module without depositing enough energy to be vetoed. 
In the fit to extract Re(e'/e) (Sec. fVC), the median of 
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the distribution of K — > irir decays inside the regenera- 
tor is used to define a perfectly sharp "effective edge," 
z e ff- For the neutral decay mode, z c g is calculated using 
the known geometry and regeneration properties, and is 
(—6.2 ± 0.1) mm from the downstream end of the re- 
generator; it is shown by the arrow above the lead in 
Fig. 0b. For the charged decay mode, z e g is determined 
by the veto threshold in the last regenerator module. The 
threshold is measured using muons collected with a sep- 
arate trigger, and results in a K — > 7r + 7r - edge that is 
(—1.65 ± 0.45) mm from the downstream edge of the re- 
generator; it is shown by the arrow above the scintillator 
in Fig. |?]b. The uncertainty comes from the geometry 
of the PMT-scintillator assembly and the threshold mea- 
surement from muons. 



4- The Veto System 

The KTeV detector uses veto elements to reduce trig- 
ger rates, to reduce backgrounds, and to define sharp 
apertures and edges that limit the detector acceptance. 
The regenerator veto was discussed in the previous sec- 
tion. 

Nine lead-scintillator (16 Xq) photon veto counters 
are positioned along the beam-line, with five located 
upstream of the vacuum-window (Fig. ||a) and four lo- 
cated downstream of the vacuum- window. These nine 
veto counters detect escaping particles that would miss 
any of the drift chambers or the Csl calorimeter. An- 
other 10 Xq photon veto is placed behind the Csl to 
detect photons that go through the beam- holes; this 
"beam-hole veto" mainly suppresses background from 
Kl — > 7T 7r 7r° decays. A scintillator bank behind 4 m 
of steel (z — 192 m in Fig. |2|) is used to veto muons, 
primarily from Kl — * ir^^v decays. 

The upstream distribution of reconstructed kaon de- 
cays is determined mainly by the "Mask Anti" (MA, 
Fig Pd), which is a 16 Xq lead-scintillator sandwich lo- 
cated at z = 123 m. The MA has two 9x9 cm 2 holes 
through which the neutral beams pass. At the down- 
stream end of the detector, the Csl crystals around the 
beam-holes are partially covered by an 8.7 Xo tungsten- 
scintillator "Collar Anti" (CA, Fig |). In addition to 
defining a sharp edge, the CA veto rejects events in which 
more than 10% of a photon's energy is lost in a beam hole. 

5. Trigger and Data Acquisition 

KTeV uses a three-level trigger system to reduce the 
total rate of stored events to approximately 2 kHz, while 
efficiently collecting K — > tttt decays. The Level 1 trigger 
(LI) has no deadtime and makes a decision every 19 ns 
(corresponding to the beam RF-structure) using fast sig- 
nals from the detector. The Level 2 trigger (L2) is based 
on more sophisticated processors and introduces a dead- 
time of 2-3/iS. When an event passes the Level 2 trigger, 



(a) Photon Veto (b) Mask Anti (MA) 




FIG. 8: (a) Transverse layout for one of five photon veto 
detectors between the regenerator and the vacuum-window 
(Fig. ^). The square aperture size varies, but is roughly 1 me- 
ter, (b) MA photon veto detector (z = 123 m). The two MA 
beam hole sizes are 9x9 cm 2 and their centers are separated 
by 20 cm. For both detectors the shaded region indicates the 
active veto area. The PMTs are at the outer edges of each 
detector. 



FIG. 9: The Collar Anti (CA) is shown as the two shaded 
"picture-frame" regions that cover the inner 60% (1.50 cm) 
of the Csl crystals surrounding the beam holes. The neutral 
beams go into the page, and the two beam-hole centers are 
separated by 30 cm. Wavelength shifting fibers transmit the 
scintillation light to PMTs at the edges of the calorimeter. 



the entire detector is read out with an average deadtime 
of 15 [is. Each event is then sent to one of twenty-four 
200-MHz SGI processors running a Level 3 (L3) software 
filter. An event passing Level 3 selection is written to 
a Digital Linear Tape for permanent storage. An inde- 
pendent set of ten 150-MHz processors is used for online 
monitoring and calibration. 

For rate reduction, the most important trigger element 
is the regenerator veto, which uses the signal from the 
downstream lead-scintillator sandwich plus signals from 
3 of the 84 scintillator modules. This veto is applied 
in Level 1 triggers to reject events from the 2 MHz of 
hadrons that interact in the regenerator. After applying 
the regenerator veto, there is still a 100 kHz rate of kaon 
decays and another 100 kHz rate of hadron interactions 
in the vacuum- window and drift chambers. Additional 
trigger requirements are used to reduce this 200 kHz rate 
by about a factor of 100 to match the bandwidth of the 
data acquisition system. 
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Separate triggers are defined for the signal and 
7r°7r° modes. Each trigger is identical for the two beams 
and for both regenerator positions. The charged, neu- 
tral, and special purpose triggers are described below. 
Each section includes a brief summary of the trigger in- 
efficiency, which is defined as the fraction of events that 
pass all analysis cuts, but fail the trigger. The ineffi- 
ciencies are measured using decays collected in separate 
minimum-bias triggers. 

Charged Decay Trigger: The Level 1 trigger for charged 
decays requires hits in the two drift chambers upstream of 
the magnet, and requires hits in the "trigger hodoscope" 
located 2 m upstream of the Csl. Each drift chamber is 
required to have at least one hit in both the x and y 
views. The trigger hodoscope consists of two 5 mm thick 
scintillation planes, each with 31 individually wrapped 
counters. There are small gaps between the counters, 
representing 1.1% of the area of each scintillation plane. 
The hodoscope counters are arranged to minimize the 
impact of these gaps; a particle traversing a gap in one 
plane cannot pass through a gap in the other plane. Each 
plane has two 14 x 14 cm 2 holes to allow the neutral 
beams to pass through without interacting. The trigger 
requires 1 or more hits in both planes, and 2 or more 
hits in at least one plane; this requirement allows for one 
of the charged particles to pass through a gap between 
the counters. The two hodoscope hits must also include 
both the upper and lower "regions," as well as the left 
and right regions. The defined hodoscope regions have 
sufficient overlap to prevent losses for K — > 7r + 7r~ events, 
except for decays in which a pion passes through a scin- 
tillator gap in the central region defined by \y\ < 7 cm. 
To reduce the trigger rate from non-7r + 7r~ kaon decays, 
signals from the veto system are used to reject events at 
the first trigger level. The muon veto is used to suppress 
Kl — > ir^n^is decays, and the photon vetos downstream 
of the vacuum-window are used to reject decays with a 
photon in the final state. 

The Level 2 trigger requires that the drift chamber hits 
in the y-view be consistent with two tracks from a com- 
mon vertex; to reduce signal loss from inefficient wires, a 
missing hit is allowed. The L3 filter selects K — ► -k + tt~ 
candidates by reconstructing two charged tracks in the 
spectrometer; to reject Kl — * 7r ± e =F ^ and Kl — > 7r + 7r~7r° 
decays, the ir + ir~ mass is required to be greater than 
450 MeV/c 2 and E/p is required to be less than 0.9 for 
both tracks. For Csl calibration and detector acceptance 
studies, 1/7 of the K e3 candidates are kept by requiring 
E/p> 0.9 for one of the tracks. 

The total charged trigger inefficiency is 0.5%, and is 
mainly from the 0.3% loss due to gaps between the scin- 
tillator counters. The drift chamber requirements result 
in a 0.1% inefficiency. Accidental effects result in a 0.06% 
loss, and there is a 0.04% loss from the trigger hardware. 
The Level 3 inefficiency is 0.09%. 

Neutral Decay Trigger: The Level 1 trigger for neu- 
tral decays is based on the total energy deposited in the 



Csl calorimeter. Using a 3100 channel analog sum, the 
threshold is 24 GeV and the 10-90% turn-on width is 
7 GeV. The photon vetos downstream of the vacuum- 
window are used in the Level 1 trigger to reduce the rate 
from Kl — * tt°tt°tt decays. The Level 2 trigger uses a 
Hardware Cluster Counter |27|] to count clusters of en- 
ergy above 1 GeV in the Csl. Four clusters are required 
for K — > 7r°7r° decays, and a separate six-cluster trigger 
is pre-scaled by five to collect Kl —> tt°tt°tt decays. The 
Level 3 filter requires that the invariant mass be greater 
than 450 MeV/c 2 for both the 7r°7r° and ttW final 
states. 

The neutral energy-sum trigger inefficiency has two 
components: a 0.6% inefficiency from early accidental 
effects and an inefficiency of 4 x 10~ 5 from the trigger 
hardware. The Level 2 trigger inefficiency of 0.4% comes 
from the Hardware Cluster Counter. The Level 3 ineffi- 
ciency is 0.01%. 

Other Triggers: In addition to the triggers used to se- 
lect specific decay modes, several special-purpose triggers 
are used for monitoring. These include (i) a "pedestal" 
trigger, which reads out the entire Csl calorimeter a few 
times per 2 s spill, (ii) a laser trigger to monitor the Csl 
(Sec. [I C 2 ), (iii) a muon trigger that requires a muon ho- 
doscope signal instead of using it in veto, (iv) a charged 
mode trigger that does not use the regenerator veto, (v) 
a trigger using only the Level 1 Csl analog sum with no 
Level 2 requirement, and (vi) an "accidental" trigger to 
record random activity in the detector that is propor- 
tional to the instantaneous intensity of the proton beam; 
triggers iii-vi are heavily pre-scaled. For the accidental 
trigger, we use a telescope consisting of three scintilla- 
tion counters, each viewed by a photomultiplier tube. It 
is located 1.8 m from the BeO target and is oriented at 
an angle of 90° with respect to the beam axis. The tar- 
get is viewed by the counters through a 6.4 x 6.4 mm 2 
hole in the stainless steel shielding around the target. A 
coincidence of signals in all three counters generates an 
accidental trigger. 

Trigger Rates: Under nominal conditions, the total 
rates passing LI, L2, and L3 are 40, 10, and 2 kHz, 
respectively; this L3 rate corresponds to approximately 
40, 000 events written to tape each minute. The dead- 
time, which is common to all triggers, is about 33% with 
roughly equal contributions from Level 2 and readout. 



D. Data Collection 

The data used in this analysis were collected in two 
distinct "E832" periods: October-December in 1996 and 
April- July in 1997. In these periods, there were about 
200 billion kaon decays between the defining collimator 
and the Csl calorimeter, of which 5 billion events were 
written to 3000 15-Gb tapes. About 5% of the events are 
used to select the K — > 7T7t sample, and the remaining 
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95% of the events are used to understand the detector. 
The neutral mode data from 1996 and 1997 are used. For 
the charged mode, the 1996 sample is not used because 
the Level 3 rejection of delayed hits in the drift chambers 
led to a 20% signal loss that is difficult to simulate. In 
1997, the Level 3 tracking was modified to avoid this loss. 
Excluding the 1996 charged mode data has a negligible 
effect on the overall statistical and systematic uncertainty 
in Re(e'/e). 

To ensure that the final data sample is of high quality, 
periods in which there are known detector problems are 
excluded. The Csl calorimeter readout suffered DPMT 
failures at a rate of about one per day, and was the 
most significant detector problem. Each DPMT failure 
was identified immediately by the online monitoring, and 
then repaired. Calibrations were frequent enough so that 
every repaired DPMT can be calibrated offline. The fi- 
nal data sample does not include periods in which there 
were dead Csl channels or dead cells in the drift cham- 
bers. Approximately 12% of the data on tape are rejected 
because of detector problems. 



E. Monte Carlo Simulation 

The Monte Carlo (MC) simulation consists of three 
main steps. The first step is kaon generation at the BeO 
target and propagation along the beamline to the decay 
point. The second step is kaon decay into an appropri- 
ate final state, and tracing of the decay products through 
the detector. The last step is to simulate the detector re- 
sponse including digitization of the detector signals. The 
simulated event format and analysis are the same as for 
the data. 

The detector geometry used in the simulation comes 
from survey measurements and data. The survey mea- 
surements are used for the transverse dimensions and z 
locations of beamline and detector elements. The trans- 
verse offsets and rotations of the drift chambers, relative 
to the Csl calorimeter, are deter mined u sing various data 

The Mask-Anti 



samples as discussed in Section II C 1 



and Collar- Anti (Sec. II C 4) aperture sizes and locations 
are determined using electrons from Kl — > ir e^v de- 
cays. 

As a result of the high flux of kaons and neutrons in 
the KTeV apparatus, there can be underlying acciden- 
tal activity in the detector that is unrelated to the kaon 
decay. After applying veto cuts, the average accidental 
energy under each Csl cluster is a few MeV, and there are 
roughly 20 extra in-time drift chamber hits. To simulate 
these effe cts, w e use data events from the accidental trig- 
ger (Sec. II C 5) to add the underlying accidental activity 
to each generated MC event. In the Csl calorimeter and 
veto system, activity from an accidental event is added 
to the MC energy deposits in a straightforward manner. 
The procedure for including accidental activity in the 
drift chamber simulation, however, is more complicated 
because an empirical model is needed to describe how an 



accidental hit can obscure a signal hit that arrives later 
on the same wire. 



1. Kaon Propagation and Decay 

The k aon energy spectrum and the relative flux of 
K° and K° states produced at the target are based on 
a parameterization ]2£| that is tuned to match KTeV 
K — > 7r + 7r~ data. The z position of each kaon decay is 
chosen based on the calculated z distribution for the ini- 
tial K° or K° state, and accounts for interference be- 
tween Kl and Ks- 

The simulation propagates the K° and K° amplitudes 
along the beamline, accounting for regeneration and scat- 
tering in the absorbers and the regenerator. Some small- 
angle kaon scatters in the absorbers (z ~ 19 m, Fig. |l|) 
can pass through the collimator system and satisfy the 
K — > 7T7T analysis requirements. The upstream colli- 
mators are modeled as perfectly absorbing, while scat- 
tering in the defining collimators and the regenerator 
are t reated using models that are tuned to data (Sec- 
tions 
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For K — > 7r + 7r~ decays, radiative corrections due to 
inner Bremsstrahlung are included [H . We do not simu- 
late the direct emission part of the radiative spectrum 
since the invariant mass cut used in the analysis es- 
sentially eliminates this component. For the K — > 7r°7r° 
mode, only the four photon final state is considered. 

GEANT |50| is used to parameterize scattering of fi- 
nal state charged particles in the vacuum-window, he- 
lium bags, drift chambers, trigger hodoscope, and the 
steel. Electrons can undergo Bremsstrahlung, photons 
can convert to e + e~ pairs, and charged pions can decay 
into a muon and a neutrino; these secondary particles are 
traced through the detector. 



2. Simulation of the Drift Chambers 

The Monte Carlo traces each charged particle through 
the drift chambers, and the hit position at each drift 
chamber plane is converted into a TDC value. The posi- 
tion resolutions measured in data are used to smear the 
hit positions, and the inverse of the x(t) map is used to 
convert the smeared hit position into a drift time. 

The simulation includes four effects that cause drift 
chamber signals to be corrupted or lost. 

1. "Wire inefficiency" results in no in- time TDC hit. 
The inefficiency is measured in 1 cm steps along 
each wire of each chamber. The average single-hit 
inefficiency is less than 1%. Since the inefficiency 
increases with distance from the wire, the measured 
inefficiency profile within the cell is used in the sim- 
ulation. 

2. A "delayed hit" results in a hit-pair with a sum- 
of-distance (SOD) that is more than 1 mm greater 
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than the nominal cell size, and there fore do es not 
satisfy the hit-pair requirement (Sec. II C 1). The 



delayed hit probability is a few percent in the re- 
gions where the neutral beams pass through the 
drift chambers, and about 1% over the rest of the 
chamber area. The effect is modeled by distribut- 
ing primary drift electrons along the track using 
a Poisson distribution with an average interval of 
340 fim, and then generating a composite signal 
at the sense wire. A delayed hit occurs when the 
signal from the nearest ionization cluster is below 
threshold, but the composite pulse from all ioniza- 
tion clusters is above threshold. The MC threshold 
is determined empirically by matching the MC de- 
layed hit probability to data; the delayed hit prob- 
ability in data is measured in 1 cm steps along each 



3. When an "in-time accidental" hit arrives before a 
signal hit on the same wire, the accidental hit is 
used instead of the signal hit because the tracking 
program considers only the first in-time hit. For 
roughly 0.7% of the hit-pairs, this effect causes a 
"low-SOD" that is more than 1 mm below the nom- 
inal cell size, and therefore does not pass the hit- 
pair requirement. An "early accidental" hit prior 
to the in-time window can also obscure a signal hit 
on the same wire for two reasons. First, the dis- 
criminator has a deadtime of 42 ns during which 
the wire is 100% inefficient. Second, large analog 
pulses can stay above the discriminator threshold 
longer than 42 ns; the variation of the pulse- length 
is modeled and tuned to data. 

4. Delta rays also cause a low-SOD for 0.5% of the hit- 
pairs. In the Monte Carlo, delta rays are generated 
in the same cell as the track, and the rate is tuned 
to match the low-SOD distribution in data. 

The quality of the drift chamber simulation is illus- 
trated by the SOD distribution in Fig. ||. Both the low 
and high-side tails in the data are well simulated. 



Simulation of the Csl Calorimeter 



the edge; this binning matches the variation in recon- 
structed position resolution, which is better for particles 
incident near the edge of a crystal. Outside the 13 x 13 
array, a GEANT-based parameterization is used instead 
of a library of individual showers; this parameterization 
models energy deposits in a 27 x 27 array. Energy leakage 
across the beam-holes is modeled based on electron data 
from Kl — > Tr^e^v decays. 

To simulate the DPMT response, the energy deposit 
in each crystal is distributed among six consecutive RF 
buckets according to the measured time profile of the 
scintillation light output. In each 19 ns wide RF bucket, 
the energy is smeared to account for photo-statistics, 
and random activity from an accidental trigger is added. 
Each channel is digitized using the inverse of the cali- 
brations obtained from data. A channel is processed if 
the digitized signal exceeds the 4 MeV readout threshold 
that was applied during data-taking. 

In addition to simulating electromagnetic showers, we 
also simulate the calorimeter response to charged pions 
and muons. The energy deposits from charged pions 
are based on a library of GEANT-based showers using 
a 50 x 50 array of small crystals. A continuous energy 
distribution is generated in x, y position bins with 4 mm 
separation. For muons, the average Csl energy deposit is 
320 MeV, and is simulated using the Bethe-Bloch energy 
loss formula. 



4- Simulation of the Trigger 

The KTeV Monte Carlo includes a simulation of the 
Level 1 and Level 2 triggers, and the Level 3 software fil- 
ter. For the K — > n + Tt~ trigger, the most important ef- 
fect to simulate is the 0.3% inefficiency due to scintillator 
gaps in the hodoscope just upstream of the Csl calorime- 
ter. The gap sizes and positions are measured in data 
using the K e % sample. The simulation also includes the 
drift chamber signals at both Level 1 and Level 2. For the 
K — > 7r°7r° trigger, we use Kl — > n ± e T i l decays to deter- 
mine the calorimeter energy-sum threshold and turn-on 
width, and to measure the Hardware Cluster Counter 
threshold for each Csl channel. 



The Monte Carlo simulation is used to predict the en- 
ergy deposit in each crystal when kaon decay products hit 
the calorimeter. In particular, the MC is needed to model 
energy leakage in the beam-holes and at the outer edges 
of the calorimeter, and to model nearby showers that 
share energy. A library of GEANT-based [H electron 
and photon showers is used to simulate electromagnetic 
showers in the calorimeter. In the shower generation, 
each electron or photon is incident on the central crystal 
of a 13 x 13 array of small crystals (32.5 x 32.5 cm 2 ). 
The showers are generated in 6 energy bins from 2 to 64 
GeV, and in x, y position bins. The position bin spac- 
ing varies from 7 mm at the crystal center to 2 mm at 



III. DATA ANALYSIS 

The analysis is designed to identify K — > 7T7t decays 
while removing poorly reconstructed events that are dif- 
ficult to simulate, and to reject background. The fol- 
lowing sections describe the analysis and the associated 
systematic uncertainties in Re(e'/e). When discussing 
systematic uncertainties, we typically estimate a poten- 
tial shift s ± cr s , where s is the shift in Re{e/e) and a s 
is the accompanying statistical uncertainty. We convert 
the shift to a symmetric systematic error, A s , such that 
the range [— A s ,+A s ] includes 68.3% of the area of a 
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Gaussian with mean s and width <j s 



a s V2iT J - 



-A. 



dx exp 



2a 2 



0.683 



(5) 



Note that A s 
s + a s /2. 



<j s when s — 0; when s > er s , A s 



A. Common Features of tt + it and 7r tt ' Analyses 

Although many details of the charged and neutral de- 
cay mode analyses are different, several features are com- 
mon to reduce systematic uncertainties. For each decay 
mode, the same cuts are applied to decays in the vac- 
uum and regenerator beams, so that most systematic 
uncertainties cancel in the single ratios used to measure 

\V+-/P\ 2 and IWH 2 - 

Since the regeneration amplitude p depends on the 
kaon momentum, we select an identical 40-160 GeV/c 
kaon momentum range for both the charged and neu- 
tral decay modes. The 40 GeV/c cut is chosen because 
of the rapidly falling detector acceptance at lower kaon 
momenta; the higher momentum cut is a compromise be- 
tween slightly higher statistics and target-i^s contamina- 
tion. We also use the same 2-vertex range of 110-158 m 
for each decay mode. The 110 m cut is chosen to be 
well upstream of the Mask Anti and therefore removes 
very few decays in the vacuum beam; this cut removes 
no events in the regenerator beam. The downstream z- 
vertex requirement avoids background from beam inter- 
actions in the vacuum- window. 

To simplify the treatment of background from kaons 
that scatter in the regenerator, the veto requirements 
for the charged and neutral mode analyses are as sim- 
ilar as possible. In both analyses, the main reduction 
in background from regenerator scattering comes from 
the requirement that there be less than 8 MeV in ev- 
ery regenerator module (Sec. II C 3 ). Interactions in the 
regenerator that are close in time to a K — > nn decay 
can add significant activity in the detector, resulting in 
events that are difficult to reconstruct. To avoid this 
problem, the trigger signal from the regenerator is used 
to reject events in a 57 ns wide window (3 RF buckets), 
which removes events immediately following or just prior 
to an interaction in the regenerator. In addition to the 
regenerator veto, we require that there be less than 150 
(300) MeV in the photon vctos upstream (downstream) 
of the vacuum-window, and less than 300 MeV in the 
Mask Anti; note that the photon veto thresholds refer to 
equivalent photon energy. 



B. K — > 7r 7r Reconstruction and Selection 

The strategy to identify K — > tt + tt~ decays is to re- 
construct two well-measured tracks in the spectrometer, 
and to reduce backgrounds with particle identification 
and kinematic requirements. 



The spectrometer reconstruction begins by finding y- 
tracks using all four drift chambers. In the x-view, sep- 
arate segments are found upstream of the magnet using 
DC1 and DC2, and downstream of the magnet using DC3 
and DC4. The extrapolated upstream and downstream 
x-track segments typically match to within 0.5 mm at the 
center of the magnet. To reduce sensitivity to multiple- 
scattering and magnetic fringe fields between the drift 
chambers, only a loose match of 6 mm is required at 
the magnet. If two x-tracks and two y-tracks are found, 
the reconstruction continues by extrapolating both sets 
of tracks upstream to define an x-z and y-z vertex. The 
difference between these two vertices, Az v t x , is used to 
define a vertex-^ 2 , 



Xltx = (AZvtx/<JAzY 



(6) 



where ctaz is the resolution of Az v t x ■ This resolution de- 
pends on momentum and opening angle, and accounts 
for multiple scattering effects. The two x-tracks and two 
y-tracks are assumed to originate from a common vertex 
if xltx < 100; this xltx requirement is sufficiently loose 
to remain insensitive to the tails in the multiple scatter- 
ing distribution. At this stage, the x and y tracks are 
independent. To determine the full particle trajectory, 
the x and y tracks are matched to each other based on 
their projections to Csl clusters; the track projections to 
clusters must match within 7 cm. 

After combining the x and y tracks, the reconstructed 
z-vertex resolution is about 30 cm near the regenerator 
and 5 cm near the vacuum- window. Each particle mo- 
mentum is determined from the track bend-angle in the 
magnet and a precise B-field map. 

An ideal event with two charged tracks results in 32 
hits in the four chambers. About 40% of the events have 
at least one missing hit or a hit-pair that does not re- 
construct a proper sum-of-distance (SOD) value. The 
tracking program can reconstruct events with many of 
these defects along the tracks. The overall single-track 
reconstruction inefficiency is measured to be 1% using 
Kl — * ir + 7r~ir° decays. 

An event is assigned to the regenerator beam if the 
regenerator x-position has the same sign as the x- 
coordinate of the kaon trajectory at the downstream face 
of the regenerator; the event is assigned to the vacuum 
beam if the signs are different. For Kl — > 7r ± e T ^ de- 
cays, this beam assignment definition cannot be used be- 
cause of the missing neutrino; instead we compare the 
x-coordinate of the decay vertex to that of the regenera- 
tor. 

The reconstruction described above is used for kaon 
decays with two charged tracks in the final state. The 
following selection criteria are specific to the K — > 7r + 7r~ 
channel, and are designed mainly to reduce background 
from semileptonic decays and kaon scattering. 

Requiring E/p < 0.85 for each track reduces the 
Kl — > 7r ± e =F j/ background by a factor of 1000. The 
Kl — > 7r ± /i =F j/ background is rejected largely by the 
muon veto in the Level 1 trigger. Since low momentum 
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muons may range out in the 4 m of steel before deposit- 
ing energy in the muon veto, we also require each track 
to have a momentum p > 8 GeV/c. To reject radia- 
tive K — > 7r + 7r~7 decays, events with an isolated electro- 
magnetic cluster above 1.0 GeV are removed if the clus- 
ter is at least 20 cm away from both extrapolated pion 
track positions at the Csl; the pion-photon separation re- 
quirement avoids removing events that have satellite clus- 
ters from hadronic interactions. To remove background 
from A — > and A — > p7r + decays, the higher momen- 
tum track is assumed to be from a proton (or antiproton) 
and the event is rejected if the proton-pion invariant mass 
is within 3.5 MeV of the known A mass. 

To provide additional background rejection of semilcp- 
tonic and K — > 7r + 7r~7 decays, the 7r + 7r _ invariant mass 
is required to be in the range 488-508 MeV/c 2 . Fig- 
ures nfla,b show the invariant mass distributions for the 
vacuum and regenerator beams after all other selection 
cuts. The shapes of the vacuum and regenerator beam 
distributions are nearly the same, and have an RMS reso- 
lution of ~ 1.6 MeV/c 2 . The low-mass tail is mainly due 
to the presence of K — > 7r + 7r~7 events in both beams. 
The tails in the vacuum beam distribution also include 
background from semileptonic decays. 
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FIG. 10: Invariant mass distributions for (a) vacuum-beam 
tv + ty~, (b) regenerator-beam tv + tt~ , (c) vacuum-beam 7r°7r°, 
(d) regenerator-beam 7r°7r°. All other analysis cuts are ap- 
plied. The vertical arrows indicate the accepted mass range. 



system (p) is projected back to a point (R re g) in the plane 
containing the downstream face of the regenerator, and 
Pt is the momentum component that is transverse to 
the line connecting R reg to the target. This definition 
of pt is used for both beams, and is optimized to dis- 
tinguish between scattered and unscattered kaons in the 
regenerator beam. The data and MC p\ distributions 



in both beams are shown in Fig 
p\ < 250 MeV 2 /c 2 , is shown by the arrows 



Target 



12; the selection cut, 



decay 
vertex 




FIG. 11: Illustration of the transverse momentum, pr, for a 
kaon that scatters in the regenerator and then decays at the 
point labeled "decay vertex." The line connecting the target 
to R reg represents the kaon trajectory upstream of the regen- 
erator, p is the measured momentum of the Tv + Tr~ system, 
and 2 r eg is the distance from the target to the downstream 
end of the regenerator. 
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FIG. 12: 7T + 7r _ distributions after all other K — > tv + ty~ 
analysis cuts in (a) the vacuum beam and (b) the regenerator 
beam. The background-subtracted data are shown as dots, 
and the histogram shows the MC prediction. The selection 
cut is shown by the arrows. The dotted histogram in the 
vacuum beam shows the MC prediction without the effect s of 
wire inefficiency, deiayed hits, and accidentais (Sec. [IE 2). 



Background from kaon scattering in the regenerator is 
suppressed mainly by the regenerator veto requirement. 
To reduce this background further, a cut is made on the 
7r + 7r~ transverse momentum, which also rejects events 
from kaons that scatter in the defining collimator. As 
illustrated in Fig. O, the total momentum of the 7r + 7r~ 



Next, we describe selection requirements that define 
apertures. For K — > 7r + 7r~, unlike K — > 7r°7r° decays, the 
edges of the Mask Anti and Collar Anti vetos are not use- 
ful apertures because of the difficulty in simulating the 
interactions of charged pions. For the MA, a charged pion 
can pass through a small amount of lead-scintillator near 
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TABLE I: Systematic uncertainties in f?e(e'/e) from the 
charged mode trigger and analysis. The uncertainties in the 
3rd col umn, which also appear in the systematics summary 
(Table VIII), are the quadratic sum of contributions in the 



2nd column. 



Source of 
uncertainty 


Uncertainty on Re(e/e) 
(xl(T 4 ) 


Trigger 
LI & L2 
L3 filter 


0.20 
0.54 


0.58 


Track reconstruction 

Alignment and calibration 
Momentum scale 


0.28 
0.16 


0.32 


Selection efficiency 
Pt cut 

DC efficiency modeling 
DC resolution modeling 


0.25 
0.37 
0.15 


0.47 


Apertures 
Wire spacing 
Regenerator edge 


0.22 
0.20 


0.30 



the edge without leaving a veto signal, and multiple scat- 
tering in the MA results in a poorly reconstructed vertex. 
For the CA, the difficulty is with pions that miss the CA, 
but interact in the Csl calorimeter; backscattering from 
these pion interactions can still deposit sufficient energy 
in the CA to produce a veto signal. To avoid these aper- 
tures, we require that the track projections be away from 
the physical boundaries of the veto detectors. In a sim- 
ilar fashion, the outer Csl aperture is defined by track 
projections so that the acceptance does not depend on 
the Csl energy profile from charged pions. 

If two tracks pass within the same drift chamber cell or 
adjacent cells, a hit from one track can obscure a hit from 
the other track. Since this effect is difficult to simulate, 
we use a track-separation requirement of two offset-cells 
(Fig. |3) between the tracks in both the x and y views for 
each chamber; offset-cells are used so that tracks near 
a wire, which have the poorest position resolution, are 
not near the boundary which defines the track-separation 
cut. This track-separation requirement results in an ef- 
fective inner aperture. 



C. 



Systematic Uncertainties from tt + -k Trigger, 
Reconstruction, and Selection 



In this section, we discuss the Re(e'/e) systematic un- 
certainties from the charged mode trigger, reconstruc- 
tion, and selection, which are summarized in Table Q. 
Systematic uncertainties related to background and ac- 
ceptance are discussed in Sections III F 5 and IV A . 



Trigger: The Level 1 and Level 2 trigger require- 
ments are studied with K e % events collected in a trig- 
ger based only on Csl calorimeter energy, and contribute 
a 0.20 x 10 -4 uncertainty to Re(e'/e). The L3 uncer- 



tainty is determined using a 1% sample of charged trig- 
gers that pass LI and L2, and are saved without requiring 
L3. From this sample, 3 x 10 5 events pass all K — > 7r + 7r~ 
analysis cuts. Applying the L3 requirement to this sam- 
ple results in a 2.2tr shift in the vacuum-to-regenerator 
ratio; the same analysis on MC events results in no shift 
from L3. These shifts in the vacuum-to- regenerator ra- 
tio correspond to a (0.45 ± 0.20) x 10~ 4 shift in Re(e'/e), 
which leads to an uncertainty of 0.54 x 10~ 4 using Eq. 0. 

Trajectory and Momentum Recons truction: The effects 
of detector misalignment (Sec. [I C 1 ) arc studied by eval- 
uating the change in Re(e'/e) as the following are changed 
within their measured uncertainties: transverse chamber 
offsets and rotations, non-orthogonality between the x 
and y wire planes, and the z-locations of the drift cham- 
bers. The time-to-distance calibration is varied to change 
the average SOD value within its uncertainty. A Re(e'/e) 
uncertainty of 0.28 x 10~ 4 is assigned based on these tests. 
The kaon mass is known to 0.031 McV pi] , leading to a 
momentum-scale uncertainty of 1 x 10~^and a Re(e'/e) 
uncertainty of 0.16 x 10~ 4 . 

Selection Efficiency: In the charged mode analysis, p\ 
is the only variable for which Re(e'/e) is sensitive to the 
cut value. Figure [l2| shows the p\ data-MC comparison 
in both beams after background subtraction (Sec. Ill F ), 
and also illustrates the importance of simulating details 
of the drift chamber performance. Increasing the p\ 
cut value from 250 to 500 MeV 2 /c 2 changes Re(e'/e) by 
(— 0.23±0.05) x 10 -4 , leading to a systematic uncertainty 
of 0.25 x 10 -4 . There is no further statistically signifi- 
cant variation if the "p\ cut value is increased beyond 
500 MeV 2 /c 2 . 

The uncertainty in modeling the drift chamber effi- 
ciency is related to the effects of delayed hit s and early 
accidentals. The delayed-hit probability (Sec. HE 2j ) pre- 
dicted by the MC is compared to data in various regions 
of the chambers and for different time periods. Resid- 
ual data-MC differences do not exceed 10% of the effect, 
which corresponds to a Re(e'/e) systematic uncertainty 
of0.21xl0~ 4 . 

There is a component of the early accidental ineffi- 
ciency that may not be modeled properly because the 
total TDC range covers only about 2/3 of the relevant 
early time window that is prior to the in-time window. 
A systematic error of 0.30 x 10~ 4 is assigned based on 
the change in Re(e'/e) when accidental hits from only 
half of the early time window are used to obscure simu- 
lated in-time hits. The total uncertainty from modeling 
the effects of delayed hits and accidentals is 0.37 x 10~ 4 
on Re(e'/e). 

The modeling of the drift chamber resolutions is 
checked by comparing the widths of the SOD distribu- 
tions between data and MC; they agree to within 5% 
which corresponds to a systematic error of 0.15 x 10~ 4 
on Re(e'/e). 
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Apertures: The drift chamber track-separation cut 
(Fig. ||) depends on the wire spacing, which is known 
to 20 fiva. on average [ p2[ . The bias in Re(e'/e) from vari- 
ations in the wire spacing is (— 0.16±0.12) x 10~ 4 leading 
to an uncertainty of 0.22 x 10~ 4 . The 0.45 mm uncer- 
tainty on the effective regenerator edge (Fig. 0b) leads 
to a small uncertainty in the expected number of K $ de- 
cays, and results in a systematic error of 0.20 x 10 -4 
Raffle). 



on 



D. 



Reconstruction and Selection 



The strategy to identify K — » 7r°7r° decays is to recon- 
struct four photon clusters in the Csl calorimeter that 
are consistent with coming from two neutral pions, and 
to reduce background with kinematic cuts. Since the re- 
construction of K — > 7r°7r° and Kl — > 7r°7r°7r decays is 
almost identical, most of the discussion will be valid for 
both these decay modes. 

The calorimeter reconstruction begins by determin- 
ing the energy and position of each cluster found in the 
Level 2 trigger. The neutral decay mode analysis requires 
two additional corrections that are not relevant for the 
Kl — > ir^e^iy calibration. First, a correction is applied 
to account for energy shared among nearby clusters. The 
second correction, described in detail at the end of this 
section, accounts for a small energy scale difference be- 
tween data and Monte Carlo. 

Each cluster is required to have a transverse energy 
distribution consistent with that of a photon. This re- 
quirement rejects decays that have significant accidental 
energy overlapping a photon cluster; for K — > 7r°7r° de- 
cays, it also rejects background from Kl — > 7r 7r°7r° de- 
cays in which two or more photons overlap at the Csl 
calorimeter. Sensitivity to modeling the transverse en- 
ergy distribution is reduced by requiring that all clus- 
ters be separated from each other by at least 7.5 cm. 
To reduce the dependence on the MC trigger simulation, 
the energy of each cluster is required to be greater than 
3 GeV. 

To reconstruct each neutral pion, we assume that two 
photon clusters originate from a ir° decay. In the small 
angle approximation, the distance between the 7r° decay 
vertex and the Csl calorimeter is given by 



r%2 



71 -^72 /^7T° ) 



(7) 



where r\i is the distance between the two photon clusters, 
E~ ll and E l2 are the two photon energies, and m^o is 
the known mass of the neutral pion. The number of 
possible photon pairings is 3 for K 
Kl — > 7T 7r 7r° decays. For each possible photon pairing 
we compute the quantity 



xio = 

3=1 



(f> — d 



(8) 



where N^o is the number of neutral pions (2 or 3), d 3 
is the distance between the Csl and the vertex of the 
j'th tt°, d avg is the weighted average of all the d; 1 , and 
CTj is the energy-dependent 7r°-vertex resolution, which 
is roughly 40 cm (30 cm) at the upstream (downstream) 
end of the decay region. The photon pairing with the 
minimum x^o value is used because this pairing corre- 
sponds to the best agreement of the ir° vertices. To re- 
duce the chance of choosing the wrong photon pairing, 
we require that the minimum x^o value be less than 12 
for K — > 7r°7r° decays, and less than 24 for Kl — * 7r 7r°7r 
decays. After all selection cuts, the probability that the 
minimum \ 2 o value gives the wrong pairing is 0.006% for 
K -> 7r°7r° and 0.02% for K L -> 7r 7r°7r°. The z-vertex 
used in the analysis is given by Zq s \ — d avg , where Zq & \ is 
the z-position of the mean shower depth in the Csl. For 
K — > 7r°7r° decays, the z-vertex resolution is ~ 30 cm 
near the regenerator and ~ 20 cm near the vacuum- 
window. 

The beam assignment is made by comparing the x- 
component of the "center-of-energy" to the regenerator 
x position. The center-of-energy is defined to be 



ik 



J2 E i 



(9) 



where Ei are the cluster energies, Xi and yi are the clus- 
ter coordinates at the Csl calorimeter, and the index i 
runs over the photons. The center-of-energy is the point 
at which the kaon would have intercepted the plane of 
the Csl if it had not decayed. The value of {x coc ,y coc } 
typically lies well inside the beam-holes, except for kaons 
that are scattered by a large angle in eithe r the regener- 
ator or the defining collimator (Sec. Ill F ) . The center- 
of-energy resolution is ~ 1 mm, which is much smaller 
than the beam separation. An event is assigned to the 
regenerator beam if ir coo has the same sign as the regen- 
erator x-position, or to the vacuum beam if ir coe has the 
opposite sign. 

Significant background rejection comes from the invari- 
ant mass cut, 490 < m_o_o < 
the photon veto cuts 
the decay vertex is assumed to be on the line joining 
the target and the center-of-energy, at a distance d n o 
from the calorimeter 
mass distributions 



moo < 505 MeV/c 2 , and from 
To calculate the invariant mass, 



Figures 10 



d show the invariant 
In the vacuum beam, events in the 



mass side-band regions are mostly from Kl — * tt tt tt 
decays in which two of the six photons are not detected. 
In the regenerator beam, the side-band regions include 
comparable contributions from 7r°7r° pairs produced in 
the lead of the regenerator (Section III F 1 ) , K 



events with the wrong photon pairing, and Kl — > 7r 7r 7r 
decays. 

The photon veto cuts are the same as in the charged 
mode analysis, with three additional cuts to reduce back- 
ground: (i) the energy of extra isolated EM clusters is 
required to be below 0.6 GeV instead of 1.0 GeV in the 
charged analysis, (ii) the photon-equivalent energy in the 
beam-hole veto behind the Csl must be less than 5 GeV, 
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and (iii) the photon-equivalent energy in the Collar Anti 
that surrounds the Csl beam-holes (Fig. ||) must be less 
than 1 GeV. 

Five apertures define the K — > ir°ir° acceptance. 

1. The Csl inner aperture near the beam holes is de- 
fined by the CA (Fig.|). 

2. The outer Csl aperture is defined by rejecting 
events in which a photon hits the outer-most layer 
of crystals (Fig. ||). This cut is applied based on 
the location of the "seed" crystals that have the 
maximum energy in each cluster. 

3. The upstream aperture in the vacuum beam is de- 
fined by the MA (Fig. |b). 

4. The upstream edge in the regenerator beam is de- 
fined by the lead at the downstream end of the 
regenerator (Figure fj]b). 

5. The requirement of at least 7.5 cm between photons 
results in an effective inner aperture. 

Since we do not measure photon angles, the transverse 
momentum is unknown and therefore cannot be used to 
reject events in which a kaon scatters in the collimator or 
regenerator. Instead, we use the center-of-energy (Eq. |^) 
to define a variable called "Ring Number:" 



RING = 40000 x Max(Aa4 e , Ay 2 oe ) , 



(10) 



where Ax coe (Ay coe ) is the x-distance (y-distance) of the 
center-of-energy from the center of the closest beam hole. 
The RING value is the area, in cm 2 , of the smallest square 
that is centered on the beam hole and contains the point 
{x C oe, Ucoe}- The vacuum and regenerator RING distribu- 
tions are shown in Fig. O. The beam size at the Csl 
is approximately 9x9 cur, corresponding to RING val- 
ues less than 81 cm 2 . The RING distribution between 
100 and 150 cm 2 is sensitive to scattering in the up- 
stream beryllium absorbers (z ~ 19 m, Fig |l|). We re- 
quire RING < 110 cm 2 , which is a compromise between 
background reduction and sensitivity to the beam halo 
simulation. 

The final part of the K — > tt°tt° analysis is to match 
the photon energy scale for data and Monte Carlo. This 
matching is performed after the background subtraction 
that will be described in Sec. |III F| . Since the energy scale 
affects the determination of both the kaon energy and the 
z- vertex position (through Eq. |?j), events can migrate into 
and out of the selected event sample depending on the 
energy scale. The final energy scale is adjusted to match 
the data and MC reconstructed z-vertex distribution of 
K — > 7r°7r° decays at the regenerator edge, as shown in 
Fig. 14i,b. Using ~ 10 6 events near the regenerator edge, 



the energy scale is determined in 10 GeV wide kaon en- 
ergy bins, which is the binning used to extract Re(e'/e). 
On average, the regenerator edge in data lies ~ 5 cm up- 
stream of the MC edge as seen in Fig. [L4b; the data-MC 
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FIG. 13: Vacuum and regenerator RING distributions for 
K —> 7r°7r° events after all other analysis cuts and after back- 
ground subtraction. The data (MC) are shown by the dots 
(histogram). The vertical arrows show the signal selection of 
RING < 110 cm 2 . 



regenerator edge difference varies by a few centimeters 
depending on the kaon energy. 

Since the vertex distance from the Csl is proportional 
to the cluster energy scale (Eq. [7]), the multiplicative 
energy scale correction is ~ 0.9992, which corresponds 
roughly to the —5 cm data-MC difference divided by the 
61 meter distance between the regenerator and the Csl. 
The regenerator edge shift in each 10 GeV kaon energy 
bin (Fig. |l4|c) is converted to an energy scale correction, 
and is applied to each cluster in data. 



E. 



Systematic Uncertainties from Trigger, 
Reconstruction, and Selection 



In this section, we discuss the Re(e'/e) systematic un- 
certainties from the neutral mode trigger, reconstruction, 
and selection, which are summarized in Table p[ System- 
atic uncertainties related to background and acceptance 



are discussed in Sections IIIF5 and IV A. 



Trigger: The Level 1 Csl "energy-sum" trigger is 
studied using the large Kl — ► ^e^v sample in the 
charged trigger which does not have LI Csl requirements 
(Sec. EIC5). The Level 2 cluster-counter is studied us- 
ing a half-million Kj, — > 7r 7r°7r° decays from a trigger 
that requires only Level 1. The Level 3 systematic un- 
certainty is based on 1.3 x 10 5 K — > 7r°7r° events that 
satisfy all analysis cuts and that were accepted online 
without requiring L3. The combined Re(e'/e) systematic 
uncertainty from the trigger is 0.18 x 10~ 4 . 

Cluster Reconstruction: The understanding of the clus- 
ter reconstruction in the Csl calorimeter, and in partic- 
ular the energy scale, results in the largest systematic 
uncertainty in the Re(e'/e) analysis. After matching the 
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FIG. 14: (a) Reconstructed K — > 7r f V° z-vertex near the 
regenerator edge for data (dots) using the K e z electron cali- 
bration, and for MC (histogram), (b) Same as (a) after the 
final energy scale adjustment is applied, (c) Data-MC differ- 
ence in the reconstructed z-vertex (meters) in 10 GeV kaon 
energy bins, before the final energy scale adjustment. 



data and MC at the regenerator edge, we use a vari- 
ety of other modes that have one or more 7r°s in the 
final state to check how well the data and MC match at 
other z-locations. All of these crosscheck samples are col- 
lected at the same time as K — > 7T7t decays, and therefore 
the detector conditions and calibration are precisely the 
same as for the K — > 7r°7r° sample. Data-MC compar- 
isons of the reconstructed neutral vertex, relative to ei- 
ther a charged mode vertex or the known vacuum-window 
location, are made in the vacuum beam for the following: 



e T e 7: 



the 



1. K — * vr°7r^,, where 7r^, refers to tt° 

reconstructed ir° — ► 77 vertex is compared to the 
e + e~7 vertex; 



2. Kl — > vr + 7T 7T° in which the reconstructed 7r° — > 
77 vertex is compared to the tt + tt~ vertex; 

3. 77 — > 7r°7r 7r°, where 77 mesons are produced by 
beam interactions in the vacuum-window at z = 
159 m; 

4. ir°ir° pairs produced by beam interactions in the 
vacuum- window at z — 159 m. 



TABLE II: Systematic uncertainties in Re(e'/e) from the neu- 
tral mode trigger and analysis. The uncertainties in the 3rd 
col umn, which also appear in the systematics summary (Ta- 
ble VIII), are the quadratic sum of contributions in the 2nd 



column. 



Source of 
uncertainty 


Uncertainty on Re(e'/e) 
(xHT 4 ) 


Trigger 

LI trigger 
L2 trigger 
L3 filter 


0.10 
0.13 
0.08 


0.18 


Cluster reconstruction 
Energy scale 
Non-linearity 
Position reconstruction 


1.27 
0.66 
0.35 


1.47 


Selection efficiency 
RING cut 

xlo cut 

Transverse shape 


0.24 
0.20 
0.20 


0.37 


Apertures 
Collar Anti 
Mask Anti 
Reg edge 
Csl size 

Photon separation 


0.42 
0.18 
0.04 
0.15 


0.48 



Additional crosschecks at the regenerator edge include 
7r 7r° pairs produced by neutron interactions in the regen- 
erator, and K* — > ir°Ks with K$ — > 7r + 7r~ . A summary 
of the neutral vertex crosschecks is given in Fig. 15 a; they 



are all consistent with the nominal energy scale correc- 
tion except for the vacuum- window 7r°7r° pairs. 

Since the vacuum-window 7r°7r° pairs result in the 
largest discrepancy, a discussion of the analysis of these 
events is given here. The selection of vacuum-window 



7r°7T pairs is similar to that for the K 



7r°7r° decay 



mode; the main difference is that K — > 7r u 7r u decays are 
excluded by selecting events in which the 7r°7r° invariant 
mass is at least 15 MeV/c 2 away from the kaon mass. The 
vacuum-window ir°ir° sample consists of 45000 events 
in the vacuum beam. The vacuum-window z-location 
is known with 1 mm precision using charged two-track 
events produced in the vacuum- window. The simulation 
of 7r°7r° pairs is tuned to match data and MC distri- 
butions of reconstructed energy, RING, and 7r°7r° invari- 
ant mass. Figure [l5| b shows the reconstructed vacuum- 
window 7T°7r vertex for data and MC after the re- 
generator edge z-distribution has been matched. The 
data-MC comparison is complicated by the helium and 
drift chamber immediately downstream of the vacuum- 
window since this extra material is also a source of 7T 7r 
pairs. The production of 7r°7r° pairs is simulated sepa- 
rately in the vacuum-window, helium, and drift cham- 
ber. To evaluate the data-MC discrepancy in the z- 
vertex distribution, a fit is used to determine the rel- 
ative contribution from each material along with the 
data-MC vertex shift. The result of this analysis is that 
the data z-vertex distribution at the vacuum- window lies 
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[2.46±0.16 (stat)±0.33 (syst)] cm downstream of the MC 
distribution for the 1997 sample; for the 1996 sample, the 
corresponding shift is [1.89±0.24 (stat)±0.33 (syst)] cm. 
The data-MC shifts are evaluated separately in the 1996 
and 1997 samples because of the different Csl PMT signal 
integration times. 
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FIG. 15: (a) The difference between the reconstructed data 
and MC neutral vertex at different locations within the decay 
region (1997 sample). The K — » 7r°7r° mode at z — 125 m 
(solid dot) has a data-MC difference of zero by definition since 
the regenerator edge has been matched. The vacuum-window 
crosscheck at z = 159 m (solid dot) shows the largest 
discrepancy; the open circles and Kl — > tt + it~tt band show 
other crosschecks. The hatched region shows the range of 
discrepancies covered by the assigned systematic uncertainty, 
(b) The reconstructed 7r°7r° vertex distribution relative to the 
vacuum-window for data after matching the regenerator edge 
(dots), for the default MC (histogram), and for the MC shifted 
2.46 cm downstream (dashed histogram). 



The data-MC "z-shift" at the vacuum-window is trans- 
lated into a Re(e'/e) uncertainty by introducing a Csl en- 
ergy scale distortion to data such that the data and MC 
z-vertex distributions match at both the regenerator and 
vacuum- window edges; this distortion changes Re(e'/e) 
by -1.08 x lfr 4 and -1.37 x 1CT 4 for the 1996 and 
1997 data samples, respectively. The combined system- 
atic uncertainty on Re(e'/e) is 1.27 x 10~ 4 . The energy 
scale distortion leading to this uncertainty varies linearly 



with decay vertex in both beams, and corresponds to the 
hatched region in Fig. [l^; non-linear energy scale varia- 
tions as a function of decay vertex are ruled out because 
they introduce data-MC discrepancies in other distribu- 
tions. 

Some reconstructed quantities in the analysis do not 
depend on the Csl energy scale, but are sensitive to en- 
ergy non-linearities. For example, the to^o^o peak varies 
by 0.2 MeV/c 2 for kaon energies between 40 and 160 
GeV, and the data-MC z-difference at the regenerator 
edge (Fig. |l4|c) varies by 4 cm over the same kaon energy 
range. Such data-MC discrepancies can be reproduced 
with cluster-energy distortions based on energy, angle, 
position, and pedestal shift. These distortions are not 
used in the final result, but are used to determine a sys- 
tematic uncertainty of 0.66 x 10~ 4 on Re(e'/e). 

Cluster position reconstruction is studied using elec- 
trons from K e 3 decays and comparing the position re- 
constructed from the Csl to the position of the electron 
track projected to the mean shower depth in the Csl. The 
position differences are parameterized and simulated; the 
maximum Re(e'/e) shift of 0.35 x 10~ 4 is assigned as a 
systematic uncertainty. 

Selection Efficiency: There are three variables for 
which Re(e'/e) is sensitive to the cut value, (i) Vary- 
ing the RING cut between 100 and 150 cm 2 changes 
Re(e'/e) by 0.24 x 10~ 4 (Fig. [If), (ii) Relaxing the xlo 
requirement such that the inefficiency from this cut is re- 
duced by a factor of 3 changes Re(e'/e) by 0.20 x 10~ 4 . 
(iii) Removing the transverse energy distribution require- 
ment changes Re(e'/e) by 0.20 x 10~ 4 . These three 
changes added in quadrature contribute an uncertainty 
of 0.37 x 10~ 4 to Re(e'/e). 

Apertures: The aperture uncertainties are mainly 
from the Collar- Anti (CA, Fig. |) and Mask-Anti (MA, 
Fig. |b). Their effective sizes and positions are measured 
with 100 /im precision using Kl — ► n e^u decays, result- 
ing in Re(e'/e) uncertainties of 0.42 x 10~ 4 and 0.18 x 10~ 4 
from the CA and MA, respectively. The Csl calorimeter 
size is known to better than 1 mm from surveys, resulting 
in a 0.15 x 10~ 4 uncertainty on Re(e'/e). The 0.1 mm un- 
certainty on the effective regenerator edge (Fig. 0b) leads 
to a Re(e'/e) uncertainty of 0.04 x 10~ 4 . Varying the 
minimum allowed photon separation between 5 cm and 
20 cm results in no significant change in Re(e'/ e). The to- 
tal Re(e'/e) uncertainty from the neutral mode apertures 
is 0.48 x 10~ 4 . 



F. Background to K 



r 7r and K 



Two types of background are relevant for this analy- 
sis. First, there is "non-7T7r" background from misiden- 
tification of high branching-ratio decay modes such as 
scmileptonic K e ^ and in the charged mode, and 

Kl — > 7r 7r°7r° in the neutral mode. The second type 
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of background is from kaons that scatter in the regener- 
ator or the defining collimator and then decay into two 
pions. Kaon scattering is the same for both the charged 
and neutral decay modes; it can be largely eliminated 
in the charged mode analysis using the reconstructed to- 
tal transverse momentum of the decay products, but the 
lack of a photon trajectory measurement does not allow 
a similar reduction in the neutral mode analysis. 

Regenerator and collimator scattering affect the 
K -> 7r°7r° reconstruction in two ways. First, a kaon can 
scatter at a small angle and still reconstruct within the 
same beam; this "in-beam" background has different re- 
generation properties and acceptance compared to un- 
scattered kaons, which can cause a bias in Re(e'/e), as 
well as in kaon parameters such as A(f> and Am. Second, 
a kaon that undergoes large-angle scattering can traverse 
from one beam to the other beam and then decay, leading 
to the wrong beam assignment in the neutral mode anal- 
ysis; this "crossover background" is mostly K$ — > ir°ir°. 
Background from regenerator scattering is roughly ten 
times larger than from scattering in the defining collima- 
tor. 

All backgrounds are simulated, normalized to data, 
and then subtracted from the K — > ttt: signal samples. 
The background-subtraction procedure is described in 
the following sections, and the background-to-signal ra- 



tios are summarized in Table III 



1. Non-Trir Background 

Ch arged pion identification and kinematic cuts 
(Sec. HI B ) eliminate most of the charged mode non-7T7r 
background; only a 0.09% contribution from the semilep- 
tonic -K^3 and K e % processes remains to be subtracted in 
the vacuum beam. Both semileptonic background distri- 
butions are simulated, and then normalized to data using 
events reconstructed outside the invariant mass and p\ 
signal region. 

The only substantial non-7T7r background in the neutral 
mode is from Kt, — > 7r 7r°7r° with undetected photons or 
photons that have merged at the Csl calorimeter. This 
0.11% background is simulated and then normalized to 
data using sidebands in the 7r°7r° invariant mass distri- 
bution in the vacuum beam. The decay Kt, —* 7r°77 has 
a branching fraction of 1.7 x 10~ 6 plfl ; it contributes 
2 x 10 -5 background in the vacuum beam and is ignored. 
We also ignore S° — > An with A — > mr°, which con- 
tributes less than 10 -5 background. 

The charged and neutral decay modes both have 
misidentification background associated with hadronic 
production in the lead plate of the last regenerator mod- 
ule (Fig. Qb). In both the charged and neutral data 
samples, this "regenerator-hadron" production is easily 
isolated in the reconstructed z-vertex distribution at the 
regenerator edge. This background is almost entirely re- 
jected by the p\ cut in the charged decay mode; the 
remaining 10 -5 background is ignored. In the neutral 



decay mode, the regenerator- 7r°7r° background is 8 x 10 -5 
(2 x 10~ 5 ) in the regenerator (vacuum) beam, and is in- 
cluded in the background simulation. 



2. Collimator Scattering Background 

Scattering in the defining collimators is studied using 
K — > 7r + 7r~ and Kt, — > 7r + 7r~7r° decays in the vacuum 
beam. Figure [l^ shows the y vs. x distribution of the 
kaon trajectory projected back to the z position of the 
defining collimator for high p\ vacuum beam events that 
satisfy all other K — > tt + tt~ requirements. The square 
bands show kaons that scattered from the defining colli- 
mator edges before decaying. The events in Fig. [16] that 
lie outside the collimator scattering bands are mainly 
from semileptonic decays. To determine the number of 
collimator scatters accurately, the roughly 10% semilep- 
tonic component is subtracted. In the charged decay 
mode, the background from collimator scattering is only 
0.01%, and is small mainly because of the p\ cut; in the 
neutral decay mode, this background is about 0.1% and 
therefore requires an accurate simulation. 

The MC simulation propagates each kaon to the defin- 
ing collimator, and checks if the kaon strikes the colli- 
mator at either the upstream end or anywhere along the 
3 meter long inner surface. Kaons that hit the collimator 
are traced through the steel and allowed to scatter back 
into the beam. A kaon that scatters in the collimator 
is parameterized to be either pure Kg or pure Kt,, with 
the relative amount adjusted to match the z-vertex dis- 
tribution for the collimator scattering sample shown in 
Fig. 

The MC treatment of collimator scattering is the same 
in both the vacuum and regenerator beams. About 1 /3 of 
the collimator-scattered kaons hit the Mask Anti (MA, 
Fig. ||), and can then punch through and exit the MA 
as cither a Kt, or K$. Based on measurements from 
data, the MC includes a kaon punch-through probability 
of 60%, and a K$ to Kt, ratio of about 50 for kaons that 
exit the MA. 



3. Regenerator Scattering Background 

In the charged mode analysis, the regenerator scatter- 
ing background is 0.074% in the regenerator beam, and 
is not present in the vacuum beam. In the neutral mode 
analysis, the corresponding background levels are 1.13% 
in the regenerator beam and 0.25% in the vacuum beam. 

Regenerator scattering is more complicated than col- 
limator scattering, particularly in the time dependence 
of K — > hit decays resulting from the coherent Kt,-Ks 
mixture. Figure [l^ shows the observed proper decay 
time distributions for kaons that scatter with small p\ 
(2500-10 4 MeV 2 /c 2 ) and large p\ (> 10 5 MeV 2 /c 2 ), 
and for the unscattered K — ► 7r + 7r~ signal (j>\ < 
250 MeV 2 /c 2 ). Note that the proper time distribution 
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FIG. 16: y vs. x distribution of the kaon trajectory 
projected back to the z position of the defining collimator 
in the vacuum beam. The sample shown here is after all 
K — » tt + 7v~ analysis cuts, except that the p% cut is changed 
to P t > 1000 MeV 2 /c 2 . The square bands correspond to the 
edges of the denning collimators. 



depends strongly on the p\ value. Large-p 2 ^ scattering 
contributes mainly to crossover background, while small- 
p\ scattering contributes mostly to in-beam background. 
A detailed description of the regenerator scattering back- 
ground is needed for the neutral decay mode because the 
K — > 7r°7r sample includes events with p\ values up to 
about 3 x 10 4 MeV 2 /c 2 in the regenerator beam, and up 
to 5 x 10 5 MeV 2 /c 2 in the vacuum beam. 

A K — > 7T7T decay from a kaon that has scattered in 
the regenerator, referred to as a "regenerator scattering 
decay," is described in the MC using a function that is 
fit to acceptance-corrected K — > tt + it~ data after sub- 
tracting collimator scattering and semileptonic decays. 
The acceptance correction allows us to fit the "true" re- 
generator scattering decay distribution, so that the scat- 
tering simulation can be used to predict background for 
both charged and neutral mode decays. As described in 
Appendix [B|, the fit function depends on proper time, 
and kaon momentum. To remove non-7r + 7r~ back- 
ground from charged tracks produced in the regenerator 
lead (Fig. ^b), the fit excludes decays within 0.2 Kg life- 
times of the regenerator edge. To avoid the p\ tail from 



this difference allows for the two scattering processes to 
be distinguished. 

Using the fit function to simulate the K — > 7r°7r° scat- 
tering background, and normalizing to data events in the 
range 300 < RING < 800 cm 2 , we find that the neutral 
mode veto requirements suppress inelastic scattering by 
an additional 16% compared to the charged mode veto 
requirements. This 16% charged-neutral difference in the 
inelastic component corresponds to a 3% charged-neutral 
difference in the total regenerator scattering background. 



— p T 2 < 250 MeV7c (coherent) 
Pt 2 



2500-1 4 MeV 2 /c 2 (diffractive) 



p T 2 > 1x10 5 MeV 2 /c 2 (inelastic+diffr) 




FIG. 17: Acceptance-corrected distribution of the number of 
regenerator beam K — > 7r + 7r _ decays vs. the number of Ks 
lifetimes (t/rs), for kaon momenta between 40 and 50 GeV/c. 
The proper time t = corresponds to the downstream edge 
of the regenerator. The p\ ranges are indicated on the plot 
for each curve. The thick curve is for coherent K — * ir + ir~ 
decays. The thin solid curve corresponds to small p\ and is 
mostly from diffractive scattering. The dashed curve is for 
large p\ events, which are due to both inelastic and diffrac- 
tive scattering. The fluctuations are due to statistics. All 
K — > 7r + 7r~ cuts except p% have been applied. Backgrounds 
from collimator scattering and semileptonic decays have been 
subtracted. 



4- Summary of Backgrounds 



Figure [l8| shows the vacuum and regenerator beam 



coherent events, only decays with p\ > 2500 MeV /c pj, distributions after all other charged mode analysis 
are used in the fit. The fit momentum region is the same 
as in the signal analysis: 40 to 160 GeV/c. 

As discussed in Section II C 3, there are two processes 



that contribute to regenerator scattering. The first pro- 
cess is diffractive scattering, which is identical in the 
charged and neutral mode analyses because no energy 
is deposited in the regenerator or photon vetos. The 
second process is inelastic scattering, which is slighlty 
different in the two modes because of different photon 
veto requirements. To address this charged-neutral dif- 
ference, the fit function is based on a phenomenological 
model that has separate terms for diffractive and inelas- 
tic scattering. The p\ distribution is significantly steeper 
for diffractive scattering than for inelastic scattering, and 



requirements. MC simulations of the background pro- 
cesses are also shown. The regenerator beam background 
is mostly from regenerator scattering, and the vacuum 
beam background is mostly from semileptonic decays. 
Figure |l^ shows the neutral mode RING distribution in 
both beams, along with MC background simulations. 
Figure ^0| shows the neutral mode background-to-signal 
ratio (B/S) as a function of the K — > tt°tt° decay vertex. 
In the regenerator beam, the main background is from 
regenerator scattering. In the vacuum beam, the largest 
sources of background are collimator scattering between 
110 m and 125 m, crossover regenerator scattering be- 
tween 125 m and 140 m, and Kl — - ;| -" -" 
z > 140 m. 



it irir decays for 
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Table III summarizes the background levels for both 
decay modes. The charged mode background level is 
~ 10~ 3 in both beams; the neutral background level is 
1.2% in the regenerator beam and 0.5% in the vacuum 
beam. The background subtraction results in corrections 
to Re(e'/e) of —12.5 x 1CP 4 for the neutral decay mode 
and —0.2 x 10~ 4 for the charged decay mode. 
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FIG. 18: Data distribution after all other K — > it + tt~ 
selection cuts for (a) the vacuum beam and (b) the regener- 
ator beam. The Monte Carlo predictions for the background 
components are overlaid. Events with p\ < 250 MeV 2 /c 2 
(vertical arrow) are included in the final K — > tv + -k~ sample. 



5. Systematic Uncertainties from Background 

The uncertainties in Re(e'/e) resulting from the back- 
ground subtraction are shown in the last column of Ta- 
bic III. In the K — > 7r + 7r~ analysis, the background con- 
on Re(e'/e); this 



tributes an uncertainty of 0.20 x 10" 
uncertainty is based on changes in Re(e'/e) when back- 
ground rejection cuts are varied for the tt + tt~ invariant 
mass, E/p, and the minimum pion momentum. 

In the K — > 7r°7r° analysis, the background contributes 
an uncertainty of 1.07 x 10~ 4 on Re(e'/e), and is mostly 
from the 5% uncertainty on the background level for in- 
beam regenerator scattering. This background subtrac- 
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FIG. 19: Data RING distribution after all other K -* ttV se- 
lection cuts for (a) the vacuum beam and (b) the regenerator 
beam. The Monte Carlo predictions for the background com- 
ponents are overlaid. Events with RING < 110 cm 2 (vertical 
arrow) are included in the final K —* 7r°7r° sample. 



tion depends largely on modeling the acceptance for high 



P Z t K 



n^n events. 



To check our understanding of the detector acceptance 



at high pj,, we use n + n pairs from Kj, 



de- 



cays. Comparing the data and MC it + tt pip distribu- 
tions, we limit the data-MC difference to be less than 
0.5% per 10000 MeV 2 /c 2 . To convert this limit on the 
li p\ slope" into a potential bias on Re(e'/e), we weight 
the pip distribution in the neutral mode background sim- 
ulation by this slope; the resulting 0.4 x 10 -4 change in 
Re(e'/e) is included as a systematic uncertainty. 

Imperfections in the phenomenological parameteriza- 
tion and fitting of the K — > ir + ir~ scattering distribu- 
tion are estimated by comparing charged mode data to 
the MC simulation. The maximum data-MC difference 
in the scattering distribution is 3% of the background 
level, which corresponds to a 0.50 x 10~ 4 uncertainty for 
Re(e'/e). 

There is also an uncertainty in how the scattering dis- 
tribution, measured with K — > ir + Tr~ decays, is used to 
simulate backgrou nd for K — » 7r 7r decays. As men- 
tioned in Section 1IIF3, the observed charged-neutral 



difference of 3% in the regenerator scattering level is ac- 
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FIG. 20: K -> n°TY° background-to-signal (B/S) ratio vs. z- 
vertex in 1 meter bins. The regenerator (vacuum) beam is 
represented by solid dots (crosses). The vertical arrow shows 
the location of the downstream edge of the regenerator. 



TABLE III: The background-to-signal ratio (B/S) for each 
subtracted background component, for each decay mode, and 
for each beam. The Re(e'/e) systematic uncertainty associ- 
ated with each background component is listed in the column 
labeled "a sys t" . 
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FIG. 21: (a) z-vertex distribution for reconstructed 
K — > 7r + 7r~ decays for the vacuum beam (thick) and regen- 
erator beam (thin histogram), (b) z-vertex distributions for 
reconstructed K — > 7r°7r° decays, (c) Kaon momentum dis- 
tributions for reconstructed K — > -k + tx~ decays, (d) Kaon 
momentum distributions for reconstructed K — > 7r°7r° decays. 
All K — > 7T7T analysis cuts have been applied and background 
has been subtracted. 



Re(e'/e) by (+0.7 ±0.3) x 10 -4 , leading to an additional 
systematic uncertainty of 0.8 x 10~ 4 . 



The other neutral mode background sources (Table III) 
have a much smaller effect on the measurement than re- 
generator scattering. The Re(e'/e) uncertainty from all 
neutral mode backgrounds is 1.07 x 10~ 4 . 



counted for by a 16% reduction in the inelastic scattering 
component in the neutral mode background simulation 
(Appendix ^). If we ignore differences between diffrac- 
tive and inelastic scattering, and simply reduce the total 
scattering level by 3% in the neutral mode simulation, 
Re(e'/e) changes by +0.3 x 10~ 4 . We assign 0.3 x 10~ 4 
as the systematic uncertainty on Re(e'/e) to account for 
the uncertainty in the inelastic-to-diffractive ratio for 
K — > 7r°7r decays, and to account for possible charged- 
neutral differences in the p\ distribution. 

We have also checked the effect of variations in 
the analysis requirements on the background subtrac- 
tion. The most significant effect is from increasing the 
regenerator-veto threshold from 8 MeV to 24 MeV in 
both the K — > 7r + 7r~ kaon scattering analysis and the 
K — > 7r°7r° signal analysis. This change doubles the 
inelastic regenerator scattering background and shifts 



G. Analysis Summary 

The numbers of events after all event selection require- 
ments and background subtraction are given in Table IV . 
The measurement of Re(e'/e) is statistically limited by 
the 3.3 million Kl — > 7r°7r° decays. Figure |2l] shows the 
z-vertex and kaon momentum distributions for the four 
event samples. 



TABLE IV: K — + nn event totals after all analysis require- 
ments and background subtraction. 





Vacuum Beam 


Regenerator Beam 


K — > 7T + 7T 


11,126,243 


19, 290, 609 


K -» ^V 


3, 347, 729 


5, 555, 789 
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IV. EXTRACTING PHYSICS PARAMETERS 

To measure physics parameters with the event samples 
described in the previous section, we correct for detector 
acceptance and perform a fit to the data. The accep- 
tance correction and t he associat ed systematic error are 

The fitting program 



described in Sections IV A IV B 



used to extract physics resul ts from the event samples 
is described in Section IV C. Systematic uncertainties 



associated with fitting are discussed in Section IV D 



A. The Acceptance Correction 

A Monte Carlo simulation is used to determine the ac- 
ceptance, which is the fraction of K — > tttt decays that 
satisfy the reconstruction and event selection criteria. 
The very different Kl and K$ lifetimes cause a differ- 
ence between the average acceptance for decays in the 
two beams. Correcting for the acceptance difference in 
the momentum and z-vertex range used in this analysis, 
the measured value of Re(e'/e) shifts by ~ 85 x 10 -4 . 
About 85% of this correction is the result of detector 
geometry, which is known precisely from optical survey 
and measurements with data. The remaining part of the 
acceptance correction depends on detailed detector re- 
sponse and resolution in the simulation. Including acci- 
dental activity in the simulation results in corrections to 
Re(e'/e) that are about -0.9 x 10" 4 and -0.5 x 10~ 4 for 
the charged and neutral d ecay modes, respectively. As 
will be discussed in Section IV B , comparisons of data and 
MC z-vertex distributions allow us to estimate the sys- 
tematic uncertainty associated with the acceptance cor- 
rection. 

For a given range of kaon p and z, the acceptance is 
defined as 



N rec l N gen 



(11) 



where iVI e „ c 



(N*> e p) is the number of reconstructed (gen- 
erated) Monte Carlo events in the specified p, z range. 
The generated p and z ranges are slightly larger than the 
ranges used in the analysis to account for the effects of 
resolution. Fig. |22| shows the acceptance as a function 
of z for 70 < p < 80 GeV/c. The K —> nn MC samples 
used to calculate the detector acceptance correspond to 
4.7 times the K — ► 7r + 7r~ data sample and 10.4 times the 
K — > 7r°7r° data sample. The resulting statistical uncer- 
tainties on Re(e'/e) from the acceptance correction are 
0.41 x 10 -4 and 0.40 x 10~ 4 for the charged and neutral 
decay modes, respectively. 

As will be described in Section IV C, we extract 



Re(e'/e) using twelve 10 GeV/c bins in kaon momentum 
and a single, integrated z bin. This p binning reduces 
sensitivity to the momentum dependence of the detector 
acceptance and to our understanding of the kaon mo- 
mentum spectrum. The use of p bins also allows us to 
account for the momentum dependence of the regenera- 
tion amplitude. 



Systematic Uncertainty From The Acceptance 
Correction 



We evaluate the quality of the simulation by compar- 
ing the data and Monte Carlo z-vertex distributions in 
the vacuum beam, where the generated z distribution 
depends only on the well known Kl lifetime. Imperfec- 
tions in the understanding of detector size and efficiency 
would change the number of reconstructed events in a 
non-uniform way along the decay region, and would re- 
sult in a data-MC difference in the z-vertex distribution. 

The procedure for converting the data-MC vertex com- 
parison into a systematic uncertainty on Re(e'/e) is as 
follows. Since Re(e'/e) is measured in 10 GeV/c kaon 
momentum bins, we weight the number of MC events 
in each energy bin so that the data and MC kaon mo- 
mentum distributions agree. We then compare the data 
and the weighted MC z distributions, and fit a line to 
the data/MC ratio as a function of z. The slope of this 
line, s, is called an acceptance "z-slope." To a good 
approximation, a z-slope affects the measured value of 
Re(e'/e) as sAz/6, where Az is the difference of the 
mean z values for the vacuum and regenerator beam ver- 
tex distributions, and the factor 6 arises from converting 
a bias on the vacuum-to-regenerator ratio to a bias on 
Re(e'/e). Az = 5.6 m and Az = 7.2 m in the charged 
and neutral hit modes, respectively. Equation ^| is used 
to convert the bias on Re(e'/e) to a systematic uncer- 
tainty. The uncertainty in tl 31 , which affects the MC 
z-vertex distribution, contributes a negligible uncertainty 
of 0.034 x 10" 4 m" 1 to the z-slope. 



Figure g3| shows the data-MC z- vertex comparisons for 
the charged and neutral tttt decay modes, and for the high 
statistics Kl n e^v and Kl — > 7r°7r 7r modes. The 
z-slope in K L -> tt+tt" is (-0.70 ± 0.30) x 10~ 4 m~\ 
and leads to a systematic uncertainty of 0.79 x 10~ 4 
in Re(e'/e). This charged z-slope has a significance of 
2.3cr, and is mostly from the first 20% of the data sam- 
ple (i.e., data collected at the start of the 1997 run). 
The very small K e % z-slope is shown as a crosscheck, 
but is not used to set the systematic error because of 
the different particle types in the final state. To as- 
sign a systematic uncertainty for the neutral decay mode 
acceptance, we use K L — + n n°iT decays. This decay 
mode has the same particle type in the final state as 
K — > 7r°7r° decays, and the Kl — > 7r 7r 7r° reconstruction 
is more sensitive to the effects of nearby clusters, en- 
ergy leakage at the calorimeter edges, and low photon 
energies. Using a sample of 50 million reconstructed 
Kl — > 7r 7r°7r° decays in both data and MC, the z-slope 
is (+0.23 ± 0.19) x 10~ 4 m _1 leading to a neutral mode 
acceptance uncertainty of 0.39 x 10~ 4 on Re(e'/e). The 
K -> z-slope of (+0.60 ± 0.53) x 10~ 4 m" 1 is con- 
sistent with the z-slope in Kl — > 7r 7r 7r decays. The 
uncertainty in the z-dependence of the acceptance for 
each decay mode is included in the summary of system- 



atic uncertainties shown in Table VIII 
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FIG. 22: MC acceptance vs. z-vertex (2 meter bins) for 
/I — > 7T7T decays with momentum range 70-80 GeV/c. The 
solid (dotted) histogram refers to the charged (neutral) decay 
mode as indicated on the plot. 



C. Fitting Decay Distributions 

For pure Kl and K$ beams, the event yields (Ta- 
ble IV) and acceptance for each mode (Fig. ^) would 
be sufficient to determine Re(e'/e) from the acceptance- 
corrected double ratio. The regenerator, however, pro- 
duces a coherent K$-Kl mixture, so that a simple dou- 
ble ratio underestimates Re(e'/e) by a few x 10~ 4 . A 
proper treatment of the regenerator, as well as target-i^s 
in both beams, is included in a fitting program. In ad- 
dition to extracting Re(e'/e), this fitting program is also 
used to make measurements of the kaon parameters Am, 

ts, 4>-\ , and A<f>; the fitting procedure described below 

applies to both Re(e'/e) and these kaon parameters. Each 
fit has different conditions related to the z-binning and 
CPT assumptions, which are summarized in Table |y|. 

The fitting procedure is to minimize the \ 2 between 
background-subtracted data yields and a prediction func- 
tion, and uses the MINUIT |33| program. The prediction 
function (V) for each beam and decay mode is 



(12) 



where is the calculated number of K 



tttt decays 

in the specified p, z range, and A p , z is the detector accep- 
tance determined by a Monte Carlo simulation (Eq. \Tl\ j . 
Note that N™ includes full propagation of the kaon state 
from the target up to the decay point, as in the MC sim- 
ulation (Sec. pE lty . 

For all fits, the prediction function is computed in 
1 GeV/cp-bins and 2 meter z-bins. To evaluate the x 2 i n 
the Re(e'/e) fit, the event yields and prediction function 
are integrated in 10 GeV/c-widep-bins, and eachp-bin is 
integrated over the the full z-range from 110 m to 158 m. 
For the other kaon parameter fits, the event yields and 
prediction function are integrated over 10 GeV/c x 2 m 
p-z bins. 

To simplify the discussion that follows, the target--K"s 
component is ignored. For a pure Kl beam, the number 
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FIG. 23: Comparison of the vacuum beam z distributions for 
data (dots) and MC (histogram). The data-to-MC ratios on 
the right are fit to a line, and the z-slopes (see text) are shown. 
The neutral distributions are for the combined 1996+1997 
samples; the charged distributions are for 1997 only. 



of K — > 7T7T decays is 



A^ocJP(p)|7j| 2 e-*/^, 



(13) 



where t = vtik^z 



y)/p is the measured proper time 



relative to decays at the regenerator edge, ij = rj^ (?7oo) 

for charged (neutral) decays, T(p) is the kaon flux, and tl 
is the Kl lifetime. The vacuum beam decay distribution 
is determined by tl ; the total event yield is proportional 
to | ry | 2 and the kaon flux. 

For a pure Kl beam incident on the KTeV regenerator, 
the number of decays downstream of the regenerator is 

7V;j ex T R {p)T reg (p) [|p(p)| 2 e^ + H 2 e- t /^ + 

2|p||?7|cos(Ami + c/>p-c/> t? )e" t/r "'' 9 , (14) 

where (f> v = arg(?7), |p| and <f> p are the magnitude 
and phase of the coherent regeneration amplitude [fl2[, 



V 



avg 



0-/T3 



1/tx)/2, Fr(p) is the kaon flux up- 
stream of the regenerator, and T reg (p) is the kaon flux 
transmission through the regenerator. The prediction 
function accounts for decays inside the regenerator by 
using the effective regenerator edge (Fig. ^b) as the start 
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FIG. 24: z decay distribution of K — * tv + tt~ decays in the 
regenerator beam, for the restricted momentum range 40- 
50 GeV/c. The MC prediction (dashed) is without the in- 
terference term that is proportional to "2|p||??|" in Eq. [l4|. 



TABLE V: Fit conditions used to analyze K — > tt-k data, "z- 
binning" refers to using 2 m z-bins in the regenerator beam. 
"Assume CPT" means that Eq. ^ is a fit constraint. Free 
parameters common to all fits, but not shown in the table, 
include the regeneration parameters |/_(70 GeV/c)| and a, 
and the kaon flux in each 10 GeV/c momentum bin, T(pi) — 







Fit Conditions 


Fit 




Assume 


Free 


Type 


z binning 


CPT 


Parameters 


Re{e'/e) 


No 


Yes 


Re(e'/e) 


Am,Ts 


Yes 


Yes 


Am, rs 


0+- 


Yes 


No 


Am, ts, 4>+- 


A0 


Yes 


No 


Am, T S , 0+_ 








Re(e'/e), Im(e'/e) 



of the decay region. All three terms in Eq. [M] are impor- 
tant, as illustrated in Fig. [24|, which shows interference 
effects in the regenerator-beam z- vertex distribution. 

Next, we discuss how the various factors in Eqs. |l3|- 
|l4| are treated in the fits. The average vacuum-to- 
regenerator kaon flux ratio [T ITr) and the average re- 
generator transmission {T re a) cancel in the Re(e'/e) fit 
as explained in Appendix [A]. To account for the mo- 
mentum dependence of the regeneration amplitude in 
the fits (see below), we need to know the momentum- 
dependence of T I (TjxF reg ); it is measured from the 
vacuum-to-regenerator ratio of Kl — > Tr + ir~ir° decays, 
and is found to vary linearly by (+7.0 ± 0.7)% between 
40 and 160 GeV/c. This variation in Tl(T R T reg ) is 
mostly from the momentum dependence of the regenera- 
tor transmission, and to a lesser extent from the movable 
absorber transmission. 

The Kl lifetime is taken from The values of Am 
and Ts are fixed to our measurements (Sec. VI A[ ) for the 
Re(e'/e) fit, and are floated in the fits for the other kaon 
parameters. In the fits that assume CPT symmetry, 0_| 



and 4>oo are set equal to the superweak phase: 

<t>r, = <Psw = tan" 1 (2Am/Ar). (15) 

The final component of the prediction function is the 
regeneration amplitude. We use a model that relates p 
to the difference between the forward kaon-nucleon scat- 
tering amplitudes for K° and K°: 



P oc/_ 



; /(0)-/(0) 



P 



(16) 



where /(0) and /(0) are the forward scattering ampli- 
tudes for K° and K°, respectively, and p is the kaon 
momentum. Additional factors that contribute to p are 
described in |34[ . 

For an isoscalar target and high kaon momentum, /_ 
can be approximated by a single Regge trajectory corre- 
sponding to the uj meson. In that case, Regge theory 
predicts that the magnitude of /_ should vary with kaon 
momentum as a power law: 



|/_(p)| = |/_(70 GeV/c)| 



P 



70 GeV/c 



(17) 



The complex phase of /_ can be determined from its 
momentum dependence through an integral dispersion 
relation, with the requirement that the forward scatter- 
ing amplitudes be analytic functions. This "analyticity" 
requirement yields a constant phase for a power-law mo- 
mentum dependence: 



<t>f- 



--(2 + a). 



(18) 



In practice, the kaon-nucleon interactions in carbon 
are screened due to rescattering processes. The effects of 
screening modify the momentum dependence of |/_(p)| 
as well as its phase. Screening corrections are evaluated 
using Glauber theory formalism |37| f or diffractive 
scattering, and using various models [|38| for inelastic 
scattering. The screening corrections in the prediction 
function used in the Re(e'/e) fit result in a 10% correc- 
tion to a, and a 0.36 x 10 -4 shift in Re(e'/e). 



D. Systematic Uncertainties from Fitting 

Uncertainties from the fitting procedure are summa- 
rized in Table VI and discussed below. These uncer- 



tainties are mainly related to regenerator properties, and 
contribute 0.30 x 10~ 4 to the Re(e'/e) uncertainty. 

The uncertainty on the momentum dependence of the 
regenerator transmission corresponds to a 0.19 x 10~ 4 
uncertainty on Re(e'/e). The sensitivity to target- Kg is 



checked by floating the A' /A: flux ratio (Sec. |ll E l| ) in 
the fit; this changes the target- K$ component by (2.5 ± 
1.6)% of itself, and leads to a systematic uncertainty of 
0.12 x 10~ 4 on Re(e'/e). 
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TABLE VI: Systematic uncertainties in Re(e'/e) from fitting. 





Re(e/e) Uncertainty 


Source of Uncertainty 


(xlO J 


Regenerator transmission 


0.19 


Target-iCs 


0.12 


Am and ts 


0.11 


Regenerator screening 


0.15 


(p p (analyticity) 


0.07 


1996 vs. 1997 K S /K L flux ratio 


0.05 


rL P 


0.02 


Total 


0.30 



TABLE VII: Fitting uncertainties in Am, ts, <f>-\ , and A<j). 

"csyst" refers to systematic uncertainty. 







CT i. f Of ' 
u syst 1W1 ■ 






Source of 


Am 




<t>+- 


A0 


Uncertainty 


(xl0 6 fi,/s) 


(xlO" 12 s) 


n 


n 


Regen. transmission 


10.0 


0.020 


0.07 


0.01 


Target-iCs 


1.4 


0.017 


0.13 


0.01 


Regen. screening 


3.0 


0.020 


0.75 


0.03 


4> p (analyticity) 


8.1 


0.030 


0.25 


0.00 


TL 


0.0 


0.001 


0.00 


0.00 


Total 


13.3 


0.045 


0.80 


0.03 



The dependence of Re(e'/e) on Am and rg is 



ARe(e'/e) = (+0.05 x 10~ 4 ) x 



Am - 5261 
+ ( _0.06*10-*)xl4^ 



, (19) 



-12 



where Am and rg are in units of 10 6 ?is _1 and 10" 
respectively. Each numerator in Eq. [l^ is the differ- 
ence between the true value and the KTeV measurement 



(Sec. VIA ); the denominators are the total KTeV uncer- 
tainties on Am and ts- Since the KTeV measurements 
of Am and rg are anti-correlated, the systematic uncer- 
tainty on i?e(e'/e) due to variations in these parameters 
is 0.11 x 10~ 4 . 

There are also uncertainties in Re(e'/e) associated with 
the analyticity relation and screening correction used 
to predict the regeneration phase 4> p . It has been ar- 
gued that the analyticity assumption is good to 0.35° in 
the E773 experiment |39|, which included kaon momenta 
down to 30 GeV/c. A smaller deviation from analyticity 
is expected with the 40 GeV/c minimum momentum cut 
used in this analysis; this leads to a 0.25° uncertainty 
in (f>p and a 0.07 x 10~ 4 uncertainty in Re(e'/e). Using 
different screening models in the fit leads to a Re(e'/e) 
uncertainty of 0.15 x 10~ 4 . 

The fitting program uses the same Ks/Kr, flux ratio 
for the charged and neutral decay modes. Since the 1996 
K — ► 7r + 7r~ sample is excluded, we consider the possibil- 
ity of a change in the kaon flux ratio between 1996 and 
1997. The kaon flux ratio depends only on the physi- 
cal properties of the movable absorber and regenerator. 
The density of these two elements could change between 
the two years because of a possible few degree temper- 
ature difference, leading to a systematic uncertainty of 
0.05 x 10~ 4 on Re(e'/e). 

Although fitting uncertainties in the Re(e'/e) measure- 
ment are a small part of the total uncertainty, they are 
more signifi cant in the measurements of Am, Tg, and 
(f>^ (Tabic VII ) . Uncertainties in the regenerator trans- 
mission and in the analyticity assumption contribute the 
largest uncertainties in the measurements of Am and 
rg. The uncertainty in the regenerator screening model 



contributes an uncertainty of 0.75° in the 0_| measure- 
ment. Fitting uncertainties have a negligible effect on 
the measurement of A<f> because of cancellations between 
the charged and neutral decay modes. 



V. MEASUREMENT OF Re(e'/t) 



The KTeV measurement of i?e(e'/e) uses the 
background-subtracted K — > tt + it~ and K — > 7r 7r° sam- 
ples in the vacuum and regenerator beams, the pre- 
diction for the K — > ttit acceptances usi ng th e Monte 
Carlo, and the fitting program. Section VA presents 



the Re(e'/e) result and a summary of the systematic un- 
certainties. Section VB presents several crosschecks, in- 



cluding a "reweighting" technique which does not use a 
Monte Carlo acceptance correction. 



A. The Re(e'/e) Result 

There are 48 measured quantities that enter into the 
Re(e'/e) fit: the observed numbers of K — > 7r+7r~ and 
K — > 7r°7r° decays in the vacuum and regenerator beams, 
each in twelve 10 GeV/c wide p bins. Within each mo- 
mentum bin we use the z-integrated yield from 110 m 
to 158 m. There are 27 fit parameters including 24 
kaon fluxes, two regeneration parameters, and Re(e'/e). 
Therefore, the number of degrees of freedom in the fit 
is 48 - 27 = 21. CPT symmetry is assumed (Eq. |||), 
and the values of Am and rg are from our measurements 



described in Sec. VIA 



For the combined 1996 and 1997 datasets, we obtain 



Re(e'/e) 
1/ (70 GeV/c)| 

X 2 /dof 



(20.71 ± 1.48) x 10~ 4 
1.2070 ± 0.0003 mbarns 
-0.5426 ±0.0008 
27.6/21 , 



(20) 



where the errors reflect the statistical uncertainties. In- 
cluding the systematic uncertainty, 

Re(e'/e) = [20.71 ± 1.48 (stat) ± 2.39 (syst)] x 10~ 4 
= (20.7 ±2.8) x 10~ 4 , (21) 
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TABLE VIII: Summary of systematic uncertainties in 
Re(e'/e). Uncertainties from the charged and neutral decay 
modes are presented in more detail in Tables [i] and [o| 

Re(e'/e) Uncertainty (xlCT 4 ) 
from: 



Source of uncertainty 


K -> ?T 


7T 


V 
A — > TT TT 


Trigger 


0.58 




0.18 


Csl energy, position recon 






1.47 


Track reconstruction 


0.32 






Selection efficiency 


0.47 




0.37 


Apertures 


0.30 




0.48 


Background 


0.20 




1.07 


^-dependence of acceptance 


0.79 




0.39 


MC statistics 


0.41 




0.40 


Fitting 




0.30 




TOTAL 




2.39 





to 



where the contributions 
are summarized in Table IVIIL 



the systematic uncertainty 
The systematic uncer- 



tainties from the charged and neutral decay modes con- 
tribute 1.26 x 1CT 4 and 2.00 x 10 -4 , respectively. The 
largest uncertainties are from the Csl energy reconstruc- 
tion (1.47 x 10 -4 ), neutral mode background subtrac- 
tion (1.07 x 10~ 4 ), z-dependence of the acceptance in the 
charged decay mode (0.79 x 10 -4 ), and the charged mode 
Level 3 filter (0.54 x 10~ 4 ). 



B. Re(e'/e) Crosschecks 

1. Consistency Among Data Subsets 

We have performed several crosschecks of our result by 
dividing the K — > tttt samples into subsets and checking 
the consistency of Re(e'/e) and other parameters among 
the different subsets. Figure |25| shows the Re(e'/e) result 
in roughly month-long time periods, in each regenerator 
position, and for the two magnet polarities. These com- 
parisons all show good agreement. The first data point 
labeled 96/97a corresponds to the current analysis ap- 
plied to the sample used in our previous publication P0| ; 
this reanalysis is discussed in Appendix |c]. 

To check the dependence on kaon momentum, 12 sepa- 
rate fits are done in 10 GeV/c momentum bins. The free 
parameters in each fit are Re(e'/e) and |/-(70 GeV/c)|, 
which is proportional to the regeneration amplitude at 
70 GeV/c. The momentum dependence of the regenera- 
tion amplitude within each 10 GeV/c bin is described by 
the power-law, p ~ p a , where a is fixed to the value found 
in the nominal fit (Eq. ^). The x 2 per degree of free- 
dom is 20.7/11 for Re(e'/e) vs. p (Fig.ffja) and 5.2/11 for 
|/_ (70 GeV/c)| vs. p (Fig. ||b); the combined x 2 /dof is 
25.9/22 . The scatter of Re(e'/e) in the higher momentum 
bins is not present in the kaon parameter measurements 
(Fig. |28|), and a linear fit to Re(e'/e) vs. p has no sig- 
nificant slope (< 0.8a). The 140-150 GeV/c bin, which 
accounts for 6.3 of the \ 2 m Fig. |2fj|a, also contributes 6.7 
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FIG. 25: Re(e'/e) consistency checks. The first five data 
points correspond roughly to one-month time periods. Reg- 
Left and Reg-Right correspond to the two regenerator posi- 
tions. Magnet-I- and Magnet— correspond to the two mag- 
net polarities, with the K —* sample common to both 
results. The uncertainties shown within each category are 
independent statistical errors. 



to the x 2 m the nominal fit (Eq. ^), with nearly equal 
contributions from all four K — > tttt samples. 

Another crosscheck is that a, which describes the re- 
generation power law (Eq. [l7]), should be the same for 
both the charged and neutral decays. A separate fit in 
each decay mode results in 



a+_ = -0.5421 ± 0.0009 (stat) 
a 00 = -0.54 45 ± 0.0017 (stat) 



(22) 



which agree to within 1.2er. 

The target-ifs correction is checked in a separate 
Re(e'/e) fit that uses only those events with kaon mo- 
menta below 100 GeV/c and a z-vertex farther than 
124 m from the BeO target. This sample has a negligible 
Ks component, and is therefore described by essentially 
pure Kl beams entering the decay region. Using this 
sub-sample, the change in Re(e'/e) from the nominal re- 
sult (Eq. is (+0.85 ± 0.89) x 10" 4 , where the error 
reflects the uncorrelated statistical uncertainty. 



2. Re(e'/e) From Reweighting Technique 

As a final crosscheck of our "standard" analysis, we 
also measure Re(e'/e) using a reweighting technique that 
does not depend on a Monte Carlo acceptance correction. 
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FIG. 27: (a) K —+ ir + ir~ track illumination in the DC1 y- 
view for vacuum (dots) and regenerator (histogram) beams, 
(b) The same track illumination after applying the weight 
factor to the vacuum beam. The vacuum-to-regenerator beam 
ratios are shown for each case in (c) and (d); note the different 
vertical scales for the ratios. 



FIG. 26: Results of twelve independent fits in 10 GeV/c 
wide p bins for: (a) Re(e'/e) and (b) |/_(70 GeV/c)|. Within 
each kaon momentum bin, Re(e'/e) and |/_(70 GeV/c)| are 
extracted from the same fit. Only statistical errors are shown. 
Each dashed line shows the average of all momentum bins. 



This technique is similar to that used by the NA48 exper- 
iment [^lj . The "local acceptance" of K — > irir decays in 
each p-z bin (1 GeV/c x 2 m) is nearly identical in both 
beams, with the only difference arising from the effects of 
accidental activity. In this method, a weight is applied 
to vacuum beam events such that the regenerator and 
weighted vacuum beam events have the same statistical 
sampling of decay vertex and kaon energy. With the same 
local acceptance in the two beams as a function of p and 
z, an ideal weight function eliminates differences in the 
reconstruction efficiencies and resolutions for K — * 7T7t in 
the two beams. The weight factor, which is applied event- 
by-event, is the a priori ratio of the regenerator beam and 
vacuum beam decay rates, 



W(p,z) 



dT reg /dt(p, z) 
dT vac /dt(p, z) 



(23) 



The functions dT vac /dt and dT reg /dt are similar to those 
given by Equations [l3] and [l4|, respectively, with the mod- 
ification that they are constrained to vanish upstream of 
the regenerator. 

The z distribution in each beam, without the weight 
factor, is shown in Fig. (Hla-b, and the pion track 



y-illumination at the first drift chamber is shown in 
Fig. [27|a. The differences between the vacuum and re- 
generator beam distributions are due to the average ac- 
ceptance difference coming from the different Kl and 
Ks lifetimes. The effect of W(p, z) is that the vac- 
uum beam distributions match those in the regenerator 
beam (Fig. |27|b). The main drawback to this reweighting 
method is that the statistical uncertainty is increased by 
a factor of 1.7 compared to the standard analysis, because 
of the loss of vacuum beam events through the reweight- 
ing function. In addition, the reweighting technique is 
more sensitive to the neutral energy reconstruction be- 
cause the weight factor (Eq. ^3|) depends on kaon energy. 

The event reconstruction and selection are very simi- 
lar to that of the standard analysis. In the K — > 7r°7r° 
mode, the event selection cuts are identical. The most 
significant difference in the reweighting analysis is the 
energy scale correction, where the absolute energy scale 
is corrected as opposed to the relative data-MC scale. 
The energy scale correction is derived from the differ- 
ence between the 77 and 7r + 7r~ z-vertex reconstructed 
in Kl — » 7r + 7r~7r°. 

For the K — > tt + tt~ selection, the reweighting analysis 
differs from the standard analysis by adding a center-of- 
energy "RING" cut (Eq. |l(J). This cut is the same for 
neutral and charged events, and eliminates the need to 
correct for the effect of kaon scattering in the movable 
absorber upstream of the regenerator (Fig. fiT). 
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Re(e'/e) is extracted in a fit which compares the 
background-subtracted yields to a prediction function. 
The prediction function is the same as Eqs.§J-(L4], except 
that in the vacuum beam, both the data and the predic- 
tion function include the weight factor. As in the stan- 
dard fit, there are 48 measured inputs, which correspond 
to the numbers of K — > tttt events in each 10 GeV/c 
wide kaon momentum bin, for both decay modes and 
both regenerator positions. The free parameters in the 
fit include 24 kaon fluxes, the regeneration amplitude and 
power-law slope (Eq. |l7| ), and Re(e'/e). In the reweight- 
ing fit, the two regeneration parameters also appear in 
the vacuum beam fit functional via the reweighting func- 
tion, where they are not varied. We find that the value of 
Re(e'/e) is quite insensitive to the regeneration parame- 
ters used in the reweighting function. As in the standard 
i?e(e'/e) fit, Am and rg are fixed to the values from the 
measurements in Sec. VIA. There arc a total of 27 free 



TABLE IX: Systematic Uncertainties in Re(e'/t) from the 
reweighting analysis. 



parameters and 21 degrees of freedom in the reweighting 
method fit. 

The systematic errors for the reweighting analysis are 
shown in Table IX. In the neutral decay mode, the largest 
systematic uncertainty results from the sensitivity to the 
minimum photon cluster energy. The other large source 
of systematic uncertainty results from accidental activity. 
The level of accidental activity in the detector is slightly 
different for the two beams; this affects the local accep- 
tance differently in the vacuum and regenerator beams 
and is not accounted for by the weight factor. In the 
standard analysis, accidental effects are accounted for in 
the Monte Carlo acceptance correction. 

To determine the uncorrelated statistical uncertainty 
between the standard and reweighting analyses, a large 
number of Monte Carlo samples are generated and fit 
with both methods; the uncorrelated uncertainty results 
mainly from the effective loss in statistics in the reweight- 
ing method. The uncorrelated systematic error is mainly 
from the uncertainty in the minimum cluster energy re- 
quirement in the reweighting analysis; there are also con- 
tributions from the acceptance correction in the standard 
analysis and accidental effects in the reweighting analy- 
sis. 

The reweighting and standard analyses were both ap- 
plied to the part of the 1997 data sample that was 
not used in the previous publication, corresponding to 
roughly 3/4 of the total sample. There is good agree- 
ment between the two analyses: the difference between 
the reweighting and standard Re(e'/e) results is 





Source of 


Re(e'/e) Uncertainty 


Sample 


Uncertainty 


(xKT 4 ) 


Neutral 


Backgrounds 


1.31 




Reconstruction 


2.93 




Trigger 


0.41 




Accidental Bias 


1.46 


v^iidi geti 


oacKgi O UIIUO 


n i r 




Reconstruction 


0.87 




Trigger 


1.09 




Accidental Bias 


0.81 




Regenerator location 


0.20 


Common 


Ring Number cut 


0.70 




Reweighting Parameters 


0.30 


Total 




3.98 



VI. MEASUREMENTS OF KAON 
PARAMETERS 

The regenerator beam decay distribution allows mea- 
surements of the kaon parameters t$ and Am, and CPT 

tests based on measurements of and A(f>. The main 

difference compared to the Re(e'/e) fit is that we fit 
the shape of decay distribution instead of the integrated 
yield. All of the fits discussed below use 2 meter wide 
z-bins in the regenerator beam from 124 m to 158 m. In 
the vacuum beam, one z-bin from 110 m to 158 m is used 
to determine the kaon flux in each 10 GeV/c momentum 
bin. 

A z-binned fit increases the sensitivity to migrations 
in z. To allow for such migrations near the regenerator 
edge, we include an extra "z-shift" parameter which is 
the shift in the effective regenerato r edge relative to the 
nominal value calculated in Section II C 3. In all fits, the 



charged and neutral data are consistent with no z-shift 
at the regenerator edge. 

Systematic errors for the kaon parameter measure- 
ments are evaluated in a manner similar to the Re(e'/e) 
analysis, and are summarized in Table [xj. The sensitivity 
to z migration is most pronounced in the A</> uncertainty 
related to Csl energy reconstruction. The complicated 
z dependence of the regenerator scattering background 
(Fig. ^0|) also contributes significant uncertainties to re- 
sults obtained from z-binned fits. 



Measurement of Am and ts 



A[i?e(e'/e)] = [+1.5 ± 2.1 (stat) ± 3.3 (syst)] x 10~ 4 
= (+1.5 ±3.9) x 10~ 4 , (24) 

where the errors reflect the uncorrelated uncertainty be- 
tween the two methods. 



To measure Am and ts, we fit the charged and neutral 
modes separately and then combine results according to 
the statistical and uncorrelated systematic errors. We 
assume CPT symmetry (Eq. |l5| ) by dynamically setting 
the value of 4> v equal to the superweak phase using the 
floated values of Am and ts- The fit values of t$ and 
Am are independent of the value of i?e(e'/e). 
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TABLE X: Systematic uncertainties in Am, ts, 4>+-, and 
A(j>. u a sys t' refers to systematic uncertainty. The uncertainty 
associated with fitting is common to both decay modes. The 
"combined total" uncertainty is explained in the text. 

cr S y S t for: 



Source of 


Am 


TS 






Uncertainty ( 


xl0 6 7i/ 


s) (xlO" 12 s) 


o 


o 


K - 


> 7r + 7r _ 


Analysis: 






Trigger 


0.2 


0.004 


0.10 


0.02 


Track reconstruction 


0.6 


0.032 


0.02 


0.02 


Selection efficiency 


3.2 


0.011 


0.35 


0.06 


Apertures 


2.8 


0.038 


0.76 


0.09 


Background 


0.8 


0.002 


0.01 


0.01 


Acceptance 


1.2 


0.026 


0.14 


0.06 


MC statistics 


2.6 


0.012 


0.28 


0.05 


Fitting 


13.3 


0.045 


0.80 


0.03 


K -> 7T+7T- Total 


14.3 


0.074 


1.20 


0.14 


K - 


Analysis: 






Trigger 


0.4 


0.013 




0.03 


Csl reconstruction 


8.1 


0.094 


- 


0.37 


Selection efficiency 


5.0 


0.035 




0.06 


Apertures 


2.2 


0.040 




0.14 


Background 


7.0 


0.030 




0.14 


Acceptance 


2.0 


0.030 




0.05 


MC statistics 


3.3 


0.016 




0.06 


Fitting 


13.3 


0.045 




0.03 


K -> tt'V Total 


18.3 


0.126 




0.43 


Combined Total 


14.2 


0.069 




0.50 


TABLE XI: Am and 


ts results for the rej 


Generator 


beam 



charged and neutral data samples, 
statistical; the second is systematic. 



The first uncertainty is 



decay 
mode 


Am 
(xlO 6 hs- 1 ) 


T5 

(xHT 12 s) 


X 2 /dof 


TY + TY~ 
TT 7T 


5266.7 ± 5.9 ±14.3 
5237.3 ±10.6 ±18.3 


89.650 ± 0.028 ± 0.074 
89.637 ±0.050 ±0.126 


228/199 
195/199 



For each charged and neutral mode fit, there are 216 
measured input quantities. The number of vacuum beam 
decays in each 10 GeV/c momentum bin gives 12 in- 
puts; the number of regenerator beam decays in each 
2 m x 10 GeV/c bin adds 17 x 12 = 204 inputs. The 
floated parameters include the kaon flux in each of 12 
momentum bins, the magnitude and phase of the regener- 
ation amplitude, a z-shift parameter, Am, and ts; these 
17 floated parameters lead to 199 degrees of freedom. The 
results of separate fits to the charged and neutral mode 

The difference between the 



data are shown in Table XI 



charged and neutral mode results, after accounting for 
the common systematic uncertainty described below, is 
1.6(7 for Am and O.ler for rg. 

Systematic errors arising from data analysis are larger 
in the neutral decay mode than in the charged mode, pri- 
marily because of larger background and uncertainties in 
the Csl energy reconstruction. The systematic uncer- 



tainties due to regeneration properties (screening, atten- 
uation, and analyticity) are more significant than in the 
Re(e'/e) analysis because there is no cancellation between 
charged and neutral mode data. These uncertainties in 
the regeneration properties are common to the charged 
and neutral mode fits, and are applied to the final result 
after averaging. The common systematic uncertainty is 
13.3xl0 6 hs" 1 on Am, and 0.045X10" 12 s on r s . 

We combine the charged and neutral mode results 
weighted by the statistical uncertainty and the indepen- 
dent parts of the systematic uncertainty. The results are 



Am 



(5261 ± 15) xlO 6 hs- 1 
(89.65 ±0.07) xlO" 12 i 



which correspond to a superweak phase of 
(j) SW = (43.38 ±0.10)°. 



(25) 
(26) 



(27) 



Measurement of < 



and 



The fit for 0_| is similar to the Am-rs fit. The 

main difference is that we remove the CPT assumption 

(Eq. [if]) and float (f>^ in addition to Am and ts- The 

fit is performed to K — > -k + -r- data only. Compared to 
the Am-Ts fit, we have the same number of measured in- 
puts (216) and one additional free parameter ), for 

a total of 216 — 18 = 198 degrees of freedom. 

There is a large correlation among c/>_| , Am, and rs, 

which is illustrated in Fig. ^l| (Appendix |d|) . The Am-Ts 
correlation is much stronger than in a fit using the CPT 
assumption, and therefore results in a larger statistical 
uncertainty on Am and ts- The 4>^ statistical uncer- 
tainty in our fit is 2.4 times larger than in a fit with a 
fixed value of ts- 

The 0.76° systematic uncertainty from apertures (Ta- 
ble |x|) is mainl y from the cell separation cut at the drift 
chambers (Sec. Ill C| ) . Different screen ing m odels result 
in a 0.75° uncertainty on <^ (Table VII). The value 



of depends on the regeneration phase 4> p ; a 0.25° 

uncertainty from the analyticity assumption leads to a 

0.25° error on The total systematic uncertainty on 

<t>+- is 1.20°. 

The results of the fit arc: 



Am 

TS 

X 2 /v 



Next, 



[44.12 ± 0.72 (stat) ± 1.20 (syst)]° 
(44.12 ± 1.40)° 
[5288 ±23 (stat)] xlO 6 hs- 1 
[89.58 ±0.08 (stat)] xlO 
223.6/197 . 



-12 



(28) 



we fit the deviation from the superweak phase, 
4>^ — 4>sWi which is a direct test of CPT symme- 
try. Compared to the (f>^ fit shown above, the fit for 

4>^ — 4>sw results in slightly reduced statistical and sys- 
tematic uncertainties because the value of 4>sw is com- 
puted dynamically using the floated values of Am and 
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ts (Eq. |l5|), and is less sensitive to the correlations. The 
result of this fit is 



4>+- - (l>sw = [+0.61 ± 0.62 (stat) ± 1.01 (syst)]° 
= (+0.61 ±1.19)° , 

and the x 2 is the same as for the </>_| fit (Eq. |2£ 



(29) 



C. Measurement of A(f> 

The measurement of A<f> is performed in a simultane- 
ous fit to neutral and charged mode data. The number of 
measured inputs is 432, which is simply twice the number 
used in the Am-rs fits, since both charged and neutral 
modes are used in the same fit. The floated parameters 
include the charged and neutral kaon fluxes in each of 
12 momentum bins (12 + 12 = 24), the regeneration am- 
plitude and phase, one z-shift term in charged and one 
in neutral, the real and imaginary parts of e'/e, Am, ts, 

and 4>.{ ; these 33 floated parameters lead to 399 degrees 

of freedom. Note that the fit uses Im(e'/e) instead of Acj) 
as a free parameter (Eq. |^) . 

The fit for Acj) benefits from the cancellation of uncer- 
tainties in the regenerator properties. Also, there is little 
correlation with the other kaon parameters such as Am, 
Tg, and the phase of e. Consequently, systematic un- 
certainties due to regenerator properties are small. The 
largest systematic uncertainty of 0.37° is from the Csl 
cluster reconstruction in the neutral mode analysis. 

There is a correlation between the real and imaginary 
parts of e'/e, with a correlation coefficient of —0.565. As a 
result, the statistical uncertainty for Re(e'/e) is increased 
by about 20% compared to the standard fit that sets 
im(e'/e) = 0. 

The results of the fit without CPT assumptions are 



Im(e'/e) 
Re(e'/e) 



[-22.9 ± 12.8 (stat) ± 26.2 (syst)] x 10" 
(-22.9 ±29.1) x 10~ 4 
[+22.5 ± 1.9 (stat)] x 10~ 4 
425/398 . 



In terms of Ac/), the result is 

Acj) = [+0.39 + 0.22 (stat) ±0.45 (syst)] c 
= (+0.39 + 0.50)°. 



neutral decay modes. Allowing each parameter to have a 
slope as a function of kaon momentum, the significance 
of each slope is between 0.5<r and 1.5a, consistent with 
no momentum dependence. 

To check the dependence on proper decay time, the 
regenerator beam K — > irir samples are divided into sub- 
sets with proper time less than and greater than 3 Kg 
lifetimes relative to the regenerator edge. This test is 
also sensitive to the significant background variations as 
a function of decay vertex (Fig. p0[ ). The entire vacuum 
beam samples are used to determine the kaon flux in 
each 10 GeV/c momentum bin, and the statistical un- 
certainty from the vacuum beam data is subtracted for 
these comparisons. For the sample with decays near the 
regenerator edge, 85% of the K — > irir decay rate is from 
the K$ term that is proportional to \p\ 2 in Eq. |lj; for the 
other sample, 42% of the K — » irir decay rate is from the 
Ks term. Fig. |2^ shows consistent results between the 
two proper time ranges for the kaon parameter measure- 
ments. Note that the measurement of 4>-\ , which has 

strong correlations with Am and ts , is more sensitive to 
early decay times; the measurement of A</>, which is very 
weakly correlated with Am and ts, is more sensitive to 
later decay times. 



VII. CONCLUSIONS 

In this paper, we report an improved measurement of 
direct CP violation in the decay of the neutral kaon: 



i?e(e'/e) 



[20.7+ 1.48 (stat) ± 2.39 (syst)] x 10~ 4 
(20.7+2.8) x 10" 4 . 

(31) 



This result, which supersedes reference pOf , is consis- 
tent with the measurement from the NA48 collaboration 
H |§ Re{e'/e) = (15.3 ± 2.6) x 10~ 4 . The average 
of our result and measurements from [|l8[ |l9|, |2l[ gives 
Re(e'/e) = (17.2 ± 1.8) x 10 -4 with a 13% confidence 
level. 

In addition, we report new measurements of the Kl- 
Ks mass difference and the Ks lifetime: 



(30) 



Am = (5261+ 15)xl0 6 hs' 1 
t s = (89.65 + 0.07) xl0~ 12 



(32) 



D. Kaon Parameter Crosschecks 

The K — » irir samples are divided into various subsets, 

among which we check the consistency of Am, ts, cf>^ , 

and Acj). For all four measurements, we find good agree- 
ment between five month-long time periods, the two re- 
generator positions, and the two magnet polarities. The 

consistency of Am, ts, 4>-\ , and Acf> as a function of kaon 

momentum is shown in Fig. ^8|. There is good agreement 
among the 12 momentum bins in both the charged and 



where CPT symmetry is assumed. Although these re- 
sults are consistent with individual previous measure- 
ments used in the PDG averages |3l), our results each 
differ from the PDG averages by more than two stan- 
dard deviations. The KTeV ts measurement is consistent 
with a recent NA48 measurement |yj , and both results 
are much more precise than the PDG average. 
Finally, we measure phase differences 

<i>+--<Psw =(+0.61 + 1.19)° . , 

Acj> = ^oo -4>+- =(+0.39 + 0.50)°, 
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FIG. 28: Kaon parameters as a function of kaon momentum 
and decay mode are shown for: (a) Am from tt + -k~ , (b) Am 

from 7r°7r°, (c) rs from 7T + 7r~, (d) rs from 7r°7r°, (e) <j>^ , 

and (f) A<f>. Each dashed horizontal line is the average of all 
momentum bins, and the x 2 /dof noted on each plot is for the 
consistency relative to the dashed line. 
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FIG. 29: Dots (open circles) show the kaon parameter mea- 
surements using regenerator beam data with proper time less 
(greater) than three Ks lifetimes relative to the regenerator 
edge. For (a) Am and (b) rs, the charged and neutral mode 
results are shown separately. For each pair of measurements 
(dot and open circle), the error bars reflect statistical uncer- 
tainties and are independent. Each vertical dashed line shows 
the nominal result using all of the data. 
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which are consistent with the CPT-symmetry prediction 
of zero. These phase differences are extracted from fits in 
which Am and t$ are free parameters to avoid the CPT 
assumptions used to extract the nominal values in Eq. |3^. 
The A4> result can be expressed in terms of Im(e'/e): 



Im(e'/e) = (-22.9 ± 29.1) x 10" 
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APPENDIX A: PRINCIPLES OF THE KTEV 
Re(e'/e) MEASUREMENT 

A simplified treatment of the KTeV measurement tech- 
nique is presented here to illustrate some of the important 
cancellations that reduce systematic uncertainties in the 
Re(e'/e) measurement. The measurement is based on the 
number of reconstructed K — > 7r + 7r~ and K — > 7r°7r° de- 
cays in the vacuum and regenerator beams. To illustrate 
how these four quantities are related to Re(e'/e), it is con- 
venient to ignore the Kg-Kh interference in both beams. 
In this simplified case, the four measured quantities are 
related to experimental parameters as follows: 



N(vaC 7T+7T ) 

N(reg 7r + 7r~) 
N(vac 7r°7r°) 
N(reg 7r°7r°) 



TS 



Br^R^R~\p\ T re g (A2) 



T L R 00-r-00 /lOOi |2 L 
B S -TV A V TOO 

TS "fCT L 

BfF™A™\p\ 2 Tre 9 , 



(A3) 
(A4) 
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where the A^'s are the observed number of K — » 7T7t de- 
cays in each beam, i?g~ ( ' 00 ^ is the branching fraction of 

Ks — ► 7r + 7r~ (7r°7r°), JFy ^ 00 - ) are the vacuum beam kaon 

fluxes, J-^ ^ 00 ^ arc the kaon fluxes just upstream of the 
regenerator, the A's are the acceptances determined by 
a Monte Carlo simulation, L ~ 40 m is the length of the 
useful decay region, T reg ~ 0.18 is the kaon flux trans- 
mission through the regenerator, 7 = Ek/Mk ~ 140 is 
the average kaon boost, and p ~ 0.03 is the regeneration 
amplitude for forward scatters. The factor L/^cts is not 
present because essentially all of the Ks decay within the 
decay region. The charged mode vacuum-to-regenerator 
single ratio is 



APPENDIX B: FUNCTION FOR KAON 
SCATTERING IN THE REGENERATOR 

This appendix describes the function used to model 
the K — > 7T7T decay distribution for kaons that scatter in 
the regenerator. The function depends on decay time, 
and kaon momentum, and is also used to separate 
diffractive and inelastic contributions. The functional 
form is 



6 



iVregscat cx V Aj-e ** |fte As < + T)e 



Az,i|2 



(Bl) 



7^ 



r. 



I 



1 



(A5) 



where A,-, a 



A 



S.L 



\Pi 



and 



are the 24 fit parameters, 



and similarly the neutral mode single ratio is obtained 
with -I — replaced by 00. To get a typical value of the sin- 
gle ratio, use \n/p\ ~ 0.07, t l /t s = 580, T v /Fr ^ 2.32 
due to the movable absorber (Fig. ||), Av/Ar ~ 0.8 and 
L/jctl ~ 0.02; this gives r ~ 0.6, which shows that 
KTeV is designed to collect roughly the same statistics 
in the vacuum and regenerator beams. The statistical 
precision on Re(e'/e) is limited by the number of vac- 
uum beam K — > 7r°7r° decays; with Ty* — 2 MHz and 
A£P ~ 0.05, the rate of K L -> 7r°7r° is - 2 Hz. The 5% 
acceptance used here is defined relative to all kaons and 
includes a livetime factor of 0.7; it is therefore smaller 
than the acceptance shown in Figure ^2] that is defined 
within a specific momentum and z-vertex range. 

The desired quantity, |r?_| /?7oo| 2 j is proportional to 

the experimentally measured double ratio, r H /V o- The 

factors tl/ts, T reg , \p\, and L/jctl cancel in the double 
ratio. For the kaon flux cancellation, we use the con- 
straint that the vacuum-to-regenerator kaon flux ratio is 
the same for both the charged and neutral decay modes, 



2^s,l, and t is the proper time of the 
decay. Note that each pj is an independent regenera- 



tion amplitude for scattering. The terms in Eq. Bl can 
be roughly associated with the known properties of kaon 
scattering in carbon, hydrogen and lead. In addition to 
these 24 parameters, there are two extra parameters that 
describe the momentum dependence of the phase (</>p . ) 
and p\ slope (oj) associated with diffractive scattering 
from the lead at the downstream edge of the regenera- 
tor. 



Of the 26 parameters in Eq. Bl 



8 are fixed based 
on previously measured properties of kaon scattering. An 
additional 12 parameters are used to float the momentum 
dependence of Ar rcgsca t in 10 GeV/c bins; the momentum 
dependence for scattered kaons varies by only a few per- 
cent compared to that of unscattered kaons. The total 
number of free parameters in the regenerator scattering 
function is 18 + 12 = 30. 

The cij parameters, which describe the exponential 
p\ dependence, are used to distinguish between inelas- 
tic and diffractive scattering. The term with the broad- 
est p\ distribution has a -1 = 2.4 x 10 5 MeV 2 /c 2 , and 
is identified with inelastic scattering. The other terms 



Rf = J-v+-/Fr^ = J-voo /3~roo — 2.32. (A6) have much steeper p\ distributions, with 5000 < a 1 < 

70000 MeV 2 /c 2 and are associated with diffractive scat- 
tering. Figure [30| shows the diffractive and inelastic con- 
tributions to regenerator scattering. 

After determining the shapes of the diffractive and in- 
elastic scattering distributions, the next step is the ab- 
solute normalization of the background relative to the 
coherent K — ► irir signal. In the charged analysis, we de- 
fine N*~ t to be the number of reconstructed K — > ir + ir~ 
events with p\ > 2500 MeV 2 /c 2 , and to be the 

„2 



Note that Eq. A6 requires equal live-time for the vac- 
uum and regenerator beams; the charged and neutral 
mode live-times, which are more difficult to control ex- 
perimentally, do not have to be the same. If the beam 
collimation system results in small kaon-flux differences 
between the two neutral beams, switching the regenera- 
tor position between the two beams ensures that Eq. A6 
is satisfied. The only quantities in the double ratio that 
do not cancel are the acceptances, for which we use a 
detailed Monte Carlo simulation. 

The quantities r_| and roo can be measured at differ- 
ent times provided that the following are the same for 
both measurements: (i) the kaon transmission in both 
the movable absorber and regenerator, (ii) the regenera- 
tion amplitude p (iii) the distance between the primary 
BeO target and regenerator. 

When Equations Al-^4 are modified to account for 
Ks~Kl interference, there is no algebraic solution for 
i?e(e'/e); results are extracted from a fit described in Sec- 
tion Fvc. 



number 
(Fig 



of coherent events with pj, < 250 MeV /c 
18). The scattering level is adjusted in the sim- 

ratio is the 



ulation so that the reconstructed N* at /N 
same in the data and MC 
tained after subtracting collimator scatters and semilep- 
tonic decays. 

In the neutral mode analysis, we define similar quan- 
tities in the regenerator beam: A^ s ° c ° at is the number of 
reconstructed K — » tt°tt events with 300 < RING < 
800 cm 2 (after subtracting the other background com- 
ponents) and AT?° h is the number with RING < 110 cm 2 



coh 

Note that N^.~ t in data is ob- 
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FIG. 30: Regenerator scattering p% distribution (solid curve) 
from a fit to acceptance-corrected K — > 7r + 7r _ data. The nor- 
malization corresponds to the measured 0.074% scattering- 
to-coherent background level in the charged decay mode (Ta- 
ble [til] ). The hatched (clear) region shows the diffractive 
(inelastic) component. The total (diffractive + inelastic) 
scattering-to-coherent ratio is shown for K — > n + n~ (solid) 
and K — > (dashed). The 3% charged-neutral difference 

reflects the additional 16% suppression of inelastic scattering 
in the neutral mode analysis. 



(Fig. |l^). Using the scattering-to-coherent ratio deter- 
mined with acceptance-corrected K — > 7r + 7r~ decays, the 
simulation over-predicts the N®® at /N®® h ratio in data by 
3%; this difference results from the additional veto re- 
quirements in the neutral mode analysis, and is illus- 
trated in Fig. |o[ To match the N°^./N°° h ratio in data, 
the neutral mode scattering simulation is adjusted by a 
16% reduction in the inelastic scattering level; note that 
the diffractive scattering level is not affected by the veto 
cuts, and is therefore the same in the charged and neutral 
analyses. 

If we ignore the difference between diffractive and in- 
elastic scattering, the simulated N^ at /N°® h ratio can be 
adjusted by a 3% reduction in the total scattering level. 
Compared with the nominal 16% reduction in the in- 
elastic level, a global 3% adjustment gives equally good 
data-MC agreement in the RING distribution from 300 
to 800 cm 2 because the diffractive and inelastic distri- 
butions are very similar for pf, > 30000 MeV 2 /c 2 . In 
the coherent signal region, however, the diffractive and 
inelastic scattering distributions are different, leading to 
a 0.01% difference in the background prediction between 
these two ways of adjusting of the N^ at /N°® h ra ti°- This 
difference is used in evaluating the systematic uncertainty 
on Re(e'/e) in Section |III F 5|. 



p0[ . The reanalysis of that data sample gives 

Re(e'/e) = [23.2 ± 3.0 (stat) ± 2.9 (syst)] x 10~ 4 

= [23.2 ±4.4 ] x 10~ 4 (CI) 
(Reference [20] sample) , 

with x 2 /dof = 18.7/21. The Re(e'je) shift is -4.8 x 10" 4 
relative to the previous result in p0| , and is larger than 
the previous systematic uncertainty of 2.8 x 10~ 4 . The 
shifts in Re(e'/e) resulti ng fr om changes in the analysis 
are summarized in Table XII . The change in the regener- 
ator scattering background for K — > ir°ir° is mainly from 
correcting an error in the kaon scattering function. The 
other shifts are due to improvements and are consistent 
with the systematic u ncer tainties assigned in ppf . All of 
the changes in Table XII are from independent sources 
which arc described below. 



TABLE XII: Sources of the Re(e'/e) shift between reference 
po| and re-analysis of the same data sample. The systematic 
cr's (crgy S t) are from |20l. 



Source of 


Re(e'/e) shift 


Number 


change 


(xlO" 4 ) 


Of (Jsyst 


Regenerator scattering 


-1.7 


2.1 


Collimator scattering 


-0.2 


0.6 


Screening 


-0.3 


1.5 


Regenerator transmission 


-0.3 


1.5 


Regenerator edge 


-0.4 


1.6 


Neutral analysis 


+0.1 


0.1 


Mask Anti MC 


+0.3 


1.3 


Absorber scattering MC 


-0.6 


0.6 


KTeV Am and r s 


-0.5 


3.1° 


MC fluctuation 


-1.1 


1.0 


Total 


-4.8 


1.7 



"Reference |20| used PDG98 values and errors for Am and Tg. 



The regenerator scattering fun ction includes two 
extra scattering terms (Eq. Bl) and two extra pa- 
rameters to simulate subtle momentum-dependent 
features. 



Collimator scattering (Sec. [II F 2 ) is measured af- 
ter subtracting semileptonic decays, and the simu- 
lation includes two extra parameters to better de- 
scribe the Kl p^-dependence. 



• In the Re(e'/e) fit (Sec. IV C), a nuclear screening 
correction is used. 

• The energy dependence of the regenerator trans- 
mission is measured with four times more 



APPENDIX C: DISCUSSION OF THE PREVIOUS 
KTEV MEASUREMENT OF Re(e'/e) 



The i?e(e'/e) result reported in Section VA includes a 
full reanalysis of our previously published data sample 



• The calculation of the effective regenerator edge in 
the neutral decay mode (Fig. 0b) includes kaon de- 
cays upstream of the last regenerator-lead piece in- 
stead of considering only decays in the last piece 
of lead. The effective regenerator edge in charged 
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mode is based on a better technique to measure the 
veto threshold. 

• Additional neutral mode veto cuts on hits in the 
trigger hodoscope and drift chambers reduce sensi- 
tivity to the transverse energy cut on Csl clusters. 

• The Mask Anti simulation allows kaons and pho- 
tons to punch through the lead-scintillator outside 
the beam- hole regions (Fig. ^). 

• The absorber scattering simulation is improved to 
match beam shapes in a special run without the 
regenerator, but with the movable absorber still in 
place. 

• We use our measurements of Am and that each 
differ by more than 2a from the PDG averages . 

The changes listed above account for —3.7 x 10~ 4 of the 
shift in Re(e'/e). We attribute the remaining shift of 
— 1.1 x 10 -4 to a ltr fluctuation between the two inde- 
pendent MC samples used for the acceptance correction. 

The largest systematic uncertainty in [^0| was based 
on the z-dependence of the acceptance for Kl — > 7r + 7r~ 
decays. Unfortunately, this data-MC discrepancy still 
remains in the reanalysis. 



APPENDIX D: CORRELATIONS AMONG KAON 
PARAMETER MEASUREMENTS 

The determination of correlations among different kaon 
parameters, including correlations arising from system- 
atic uncertainties, is based on the following \ 2 '- 



The total uncertainties and correlation coefficients for 



X 2 (A,a. ( ) 



X 



i 



X + J2i Pi a i ■ 



(Dl) 



X is the vector of the measured physics quantities from 
the fit (e.g., Am and t<j), and C is the covariance ma- 
trix including statistical uncertainties for both data and 
MC. The «i are fit parameters representing the number 
of "tr" associated with each source of systematic uncer- 
tainty that causes correlations among the physics quanti- 
ties (a; = 1 corresponds to la change). Xq is the vector 
of physics quantities obtained in the nominal fit with all 
cti = 0, and Pi are the changes in the central value of X 
corresponding to a lcr change in the systematic source i. 
The total uncertai nty in a physics quantity obtained 



from the x 2 m Eq. Dl is equivalent to adding statis- 
tical and systematic uncertainties in quadrature. The 
minimization of this % 2 accounts for correlations among 
the physics quantities for each source of systematic un- 
certainty. We ignore correlations among the differ- 
ent sources of systematic uncertainty; these correlations 
are reduced by the grouping of systematic uncertainties 
shown in Table [x| 



all of the fits are given in Table XIII . Note that the sta- 
tistical and systematic uncertainties are larger in the fits 
in which CPT symmetry is not assumed. Figure |3l| shows 
la contours of statistical and total uncertainties for the 

measurements of Am, ts, and 0_| For the Am-rs fit 

(Fig. |3l|a) , systematic uncertainties have a significant ef- 
fect on the Am-rs correlation. For the fit (Fig. |3l|b- 

d), systematic uncertainties have a very small effect on 
the correlations. Figure |3^ shows the correlation between 
the real and imaginary parts of e'/e from the fit without 
CPT assumptions; systematic uncertainties have a small 
effect on the correlation. 



TABLE XIII: Total uncertainties and correlation coefficients 
for the three fits. The fit parameters and CPT assumptions 
are listed above each set of values. 



I. Am-rs Fit with CPT Assumption: 

Am 




Total Error lSxlO far 1 


0.07x10"" s 


Correlation coefficients: 

Am 1. 
r s -0.396 


1. 


II. Am-Ts-(j>+- Fit without CPT Assumption: 

Am ts </)+- 


Total Error 42x10° far 1 


0.13xl0 _l: " s 1.40° 


Correlation coefficients: 

Am 1. 
r s -0.874 
4>+- +0.987 


1. 

-0.898 1. 


III. e'/e Fit without CPT Assumption: 

Re{e'/e) 


Jm(e'/e) 


Total Error 4.0 x 10 -4 


29.1 x 10~ 4 


Correlation coefficients: 

Re(e'/e) 1. 
Im(e'/e) -0.647 


1. 
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FIG. 32: la contours of statistical uncertainty (dashed) and 
total uncertainty (solid) for Re(e'/e) vs. Im(e'/e) from the fit 
with no CPT assumptions related to </>_] — and (poo- 



FIG. 31: lcr contours of statistical uncertainty (dashed) 
and total uncertainty (solid) for (a) Am-rs fit (combined 
charged+neutral) , and for (b)-(d) correlations from the 4>^ — 
fit. 
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